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ABSTRACT 

Lukens ,  Larry  Allan,  Ph.D.,  Purdue  University, 
March  19  71.   An  Experimental  Investigation  of  Electric 
Field  Intensity  and  Wall  Heat  Transfer  for  the  Heating 
Region  of  a  Constricted  Arc  Plasma.   Major  Professor: 
Frank  P.  Incropera. 

Electric  field  intensity  and  total  wall  heat 
transfer  measurements  have  been  performed  for  the  inlet, 
asymptotic,  and  anode  regions  of  a  laminar,  atmospheric 
cascade  arc.   In  addition  to  the  diagnostics  required 
for  these  measurements,  a  unique  radiation  gage  has 
been  developed  and  used  to  obtain  radiative  wall  heat 
flux  data.   While  the  majority  of  the  measurements  have 
been  made  for  an  argon  arc,  some  data  are  reported  for 
helium  and  nitrogen  arcs. 

Comparisons  between  the  experimental  data  and 
available  equilibrium  and  nonequilibrium  numerical 
solutions  are  presented  in  an  attempt  to  ascertain  the 
adequacy  of  single-fluid  models  to  accurately  describe 
overall  arc  energetics.   For  the  most  part,  the  equili- 
brium theory  predictions  for  the  asymptotic  region  in  an 
argon  arc  agree  with  the  experimental  results  to  within 
the  combined  experimental  and  theoretical  uncertainties. 
However,  there  is  a  small,  but  discernable  nonequilibrium 
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effect  at  both  the  lower  and  higher  arc  currents  in  argon. 
In  contrast,  the  nonequilibrium  theory  predictions  are 
in  excellent  agreement  with  the  experimental  results  over 
the  entire  range  of  asymptotic  conditions  investigated  in 
argon.   Further  comparisons  are  presented  which  indicate 
that  atmospheric  helium  arcs  are  characterized  by  a  signi- 
ficant departure  from  local  thermodynamic  equilibrium. 

On  the  basis  of  comparisons  between  experiment  and 
theory  for  the  anode  region,  it  is  shown  that  anode  heat 
balance  models  are  ill-suited  for  providing  design  es- 
timates for  the  total  heat  flux  to  the  anode  of  the  arc 
constrictor. 


CHAPTER  I   INTRODUCTION 

1.1   General  Objectives  of  This  Study 
Owing  to  interest  in  space  flight  technology,  MHD 
systems,  and  thermonuclear  energy  conversion,  considerable 
effort  has  been  devoted  to  the  study  of  high  temperature 
gas  flows.   Of  particular  interest  to  this  study  is  that 
portion  of  this  effort  which  deals  with  a  device  which  has 
long  served  as  a  versatile  tool  in  generating  high  tempera- 
ture gases  for  both  industrial  and  laboratory  use  -  the 
plasma  arc  generator. 

Designed  initially  by  Maecker  ,  the  operation  of 
the  Gerdien,  or  constricted  tube,  arc  generator  is  based 
upon  the  conversion  of  electrical  to  thermal  energy.  In 
the  most  common  configuration  (Figure  1.1),  this  conversion 
is  accomplished  by  passing  a  gas  through  an  electric  arc 
which  is  established  by  applying  a  d-C  potential  between 
the  electrodes. 

The  constricted  tube  arc  produces  a  continuous  and 
well-defined  plasma  flow  which,  in  addition  to  its  use  in 
certain  technological  applications,  is  particularly  useful 
for  the  experimental  determination  of  high  temperature 
gas  transport  properties  and  the  development  of  high 


temperature  diagnostic  techniques.   Irrespective  of  its 
use,  however,  the  constricted  tube  arc  is  distinguished 
by  extremely  high  heat  fluxes  to  the  containing  walls  and 
these  walls  are  often  destroyed  if  sufficient  cooling  is 
not  provided.   Accordingly,  an  understanding  of  the  elec- 
trical and  thermal  characteristics  of  arc  constrictors 
is  important  to  the  designer  as  well  as  to  those  with  a 
general  interest  in  arc  physics.   The  objective  of  this 
study  is  to  enhance  this  understanding  through  an  experi- 
mental investigation  of  the  electrical  and  thermal  effects 
in  an  atmospheric  cascade  arc  and  through  an  interpretation 
of  these  effects  in  terms  of  available  theoretical  studies. 
While  the  major  portion  of  the  experimental  work  has  been 
performed  for  an  argon  plasma,  some  measurements  have 
also  been  made  for  nitrogen  and  helium  arcs. 

The  regions  of  interest  for  the  cascade  arc  are 
depicted  schematically  in  Figure  1.1.   The  arc  column, 
which  is  a  composite  of  individually  water-cooled  and 
electrically  isolated  copper  discs,  is  divided  into  three 
distinct  regions:  inlet,  asymptotic,  and  anode.   The  gas 
flow  in  each  region  is  subjected  to  an  externally  applied 
electric  field,  and  the  inlet  and  asymptotic  regions  are 
generally  bounded  by  non-current-carrying  surfaces.   In 
contrast,  the  anode  surface  is  electrically  conducting. 
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Figure  1.1  Designation  of  Arc-Heating  Regions 


The  inlet  region  is  further  distinguished  by  a  significant 
axial  variation  of  the  flow  variables,  whereas  the  asymp- 
totic region  is  characteristic  of  thermally  fully-devel- 
oped flow.   As  a  matter  of  convention,  the  axial  location 
at  which  the  mean  (mass  averaged)  enthalpy  attains  90% 
of  its  fully-developed  value  is  chosen  to  mark  the  be- 
ginning of  the  asymptotic  region. 

1.2   Previous  Investigations  of 
Plasma  Tube  Flow 
Although  the  literature  devoted  to  theoretical  study 
of  the  inlet  and  asymptotic  regions  is  extensive ,  few 
experimental  studies  have  been  performed  for  purposes  of 
determining  the  overall  arc  energetics.   Of  particular 

scarcity  are  experimental  studies  pertaining  to  the  inlet 

[21 
region  of  the  water-cooled  duct.   Skifstad    performed 

electrical  potential  and  wall  heat  flux  measurements  in 

the  inlet  region  of  a  helium  arc,  and  investigators  at 

r  3  —  5 1 

the  NASA-AMES  Research  Center1    J  have  made  similar 

measurements  in  air  and  nitrogen  arcs.   However,  the  range 
of  experimental  conditions  (gas  flow  rate,  arc  current, 
and  duct  geometry)  investigated  in  these  studies  differ 
significantly  from  those  of  this  investigation.   Conse- 
quently, meaningful  comparisons  cannot  be  made  between  the 
data  of  this  study  and  the  experimental  results  reported 
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by  these  investigators.   Similarly,  Pfender,  et  al.   ' 
report  results  for  the  local  wall  heat  flux  in  the  inlet 
region  of  an  argon  cascade  arc;  however  these  measurements 
were  for  turbulent  flow  and  were  used  primarily  for  the 
support  of  anode  heat  transfer  studies. 

A  greater  number  of  experimental  investigations 
have  been  performed  for  the  asymptotic  region  of  the 
arc      .   While  these  studies  are  primarily  concerned 
with  the  use  of  an  inversion  method  to  determine  plasma 
transport  properties,  they  necessarily  include  measure- 
ment of  the  arc-column  electrical  characteristics  and 

are  therefore  germane  to  this  study*.   In  addition,  Morris, 

[14] 
et  al.     present  asymptotic  wall  heat  flux  data  and 

Emmons     and  Maecker     report  on  the  percentage  con- 
tribution of  radiation  to  the  total  energy  balance. 

In  view  of  the  excellent  reviews  given  by 
Skifstad     and,  more  recently,  by  Bower     ,  the  present 
discussion  of  theoretical  investigations  is  limited  to 
those  studies  which  are  to  be  used  for  comparison  with 
experiment  in  Chapter  IV. 

The  approximate  analytical  solution  by  Stine  and 

[181 
Watson     was  one  of  the  first  reported  theoretical 

treatments  of  the  inlet  region  of  the  cascade  arc.   This 


*  The  inversion  method  of  property  determination  will  be 
discussed  further  in  Chapter  IV. 


solution  is  based  upon  a  reduced  form  of  the  governing 
energy  equation,  obtained  by  neglecting  radiation,  radial 
convection,  and  viscous  dissipation.   In  addition,  the 
Stine-Watson  model  assumes  a  uniform  mass  flux  throughout 
the  constrictor  and  uses  linearized  relationships  between 
the  enthalpy  and  related  properties  (T,  k/C  ,  and  a).   The 

IT 

primary  weakness  of  this  model  is  its  failure  to  properly 
account  for  the  variation  of  the  transport  properties 
with  temperature  and,  particularly  when  applied  to  the 
higher  temperature  argon  and  nitrogen  arcs,  its  failure 
to  account  for  radiative  losses. 

To  remove  the  need  for  the  above  assumptions,  de- 
tailed numerical  solutions  for  the  inlet  region  have  been 

[191 
carried  out  by  Watson  and  Pegot     and  by  Bower  and 

Incropera     .   Of  the  two,  the  latter  will  be  used  for 
the  comparisons  of  this  study  since  it  gives  greater  con- 
sideration to  pertinent  thermal  effects.   Specifically, 
Bower  and  Incropera  formulate  a  rigorous  equilibrium  model 
for  the  laminar,  steady,  and  axially  symmetric  (without 
swirl)  flow  of  an  atmospheric  argon  plasma  through  the 

heating  region  of  a  constricted  tube  arc.   An  implicit 

T17] 
finite-difference  scheme     is  then  used  to  solve  the 

resulting  conservation  equations  (mass,  momentum,  and 

energy)  consistent  with  the  thermal  and  caloric  equations 


of  state,  Ohm's  law,  and  the  temperature  dependence  of  the 
required  thermodynamic  and  transport  properties.   Implicit 
in  the  final  form  of  the  equations  are  the  boundary  layer 
approximations  and  the  assumptions  of  negligible  radial 
current  density,  quasi-neutrality ,  and  local  thermodynamic 
equilibrium.   The  results  of  this  numerical  solution  are 
presented  in  the  form  of  wall  heat  transfer  and  shear 
stress  (as  predicted  from  energy  and  momentum  balances)  and 
correlations  between  the  mean  Nusselt  number  and  an  Ohmic 
heating  parameter. 

As  a  result  of  the  evolution  of  the  aforementioned 
numerical  solutions,  it  would  appear  that  the  theoretical 
treatment  of  the  inlet  and  asymptotic  regions  of  arc 

constrictors  has  reached  a  practical  limit.   However,  the 

[13] 
variable  tube  diameter  experiments  performed  by  Emmons 

and,  more  recently,  the  spectroscopic  measurements  con- 

T211  f 22] 

ducted  by  Bott     and  by  Incropera  and  Giannaris 

provide  strong  experimental  evidence  of  the  presence  of 

nonequilibrium  in  confined  plasmas.   Accordingly,  the 

validity  of  the  local  thermodynamic  equilibrium  assumption 

which  is  implicit  in  the  single-fluid  model  is  subject  to 

question.* 


A  comprehensive  discussion  of  local  thermodynamic      too] 
equilibrium  in  confined  plasmas  is  given  in  references 

and  f23' 


In  an  attempt  to  account  for  the  nonequipartition 

of  species  translational  energies  (thermal  nonequilibrium) , 

[241 
Pytte  and  Winsor     formulated  a  nonequilibrium  theory 

for  the  asymptotic  region  of  a  helium  arc.   This  attempt, 
however,  is  at  best,  semi-empirical,  for  its  solution 
relies  upon  experimentally  determined  arc-column  charact- 
eristics.  In  addition,  the  Pytte-Winsor  model  suffers 
from  use  of  the  Saha  (equilibrium)  equation  to  compute 
the  plasma  composition. 

While  accounting  for  both  thermal  nonequipartition 

[251 
and  chemical  nonequilibrium,  Clark     has  formulated  a 

rigorous  multifluid  (electrons,  ions,  and  neutral  atoms) 

model  for  the  steady,  laminar  flow  of  a  nonequilibrium 

argon  plasma  in  the  inlet  and  asymptotic  regions  of  a 

constricted  tube  arc.   After  applying  the  boundary  layer 

approximations,  the  governing  set  of  species  and  global 

conservation  equations,  along  with  the  thermally  perfect 

gas  equation  of  state  and  Ohm's  law,  are  solved  numerically 

T251 
using  finite-difference  techniques     .   Of  particular 

interest  to  this  study  are  the  predictions  for  the  asymp- 
totic electric  field  intensity,  total  wall  heat  flux,  and 
ratio  of  radiative  to  total  wall  heat  flux  which  result 
from  the  nonequilibrium  solution. 

The  need  for  extended  operating  times  for  the  most 
severely  heated  component  of  the  arc  constrictor,  the 


anode,  has  prompted  considerable  research  on  the  complex 
energy  exchange  mechanisms  occurring  at  this  surface. 
Energy  transfer  models  have  been  proposed  in  which  the 
various  contributions  to  the  total  anode  heat  flux  are 
delineated.   For  the  most  part,  these  models  are  based  on 
a  heat  balance  across  the  anode  surface   '  '       *.   In 
addition,  extensive  studies  of  the  electrical  characteris- 
tics of  anodes  have  been  conducted  by  researchers  at 

[31-331 
Dartmouth       .   For  purposes  of  comparison  with  the 

anode  heat  transfer  measurements  of  this  study,  primary 

emphasis  will  be  placed  upon  the  widely  used  model  which 

r  6 1 

is  described  by  Pfender,  et  al.    . 

1.3   Scope  and  Specific  Objectives  of  This  Study 

The  primary  objective  of  this  study  is  the  measure- 
ment of  the  axial  distribution  of  the  electric  field 
intensity  and  the  total  wall  heat  flux  for  the  inlet, 
asymptotic,  and  anode  regions  of  an  atmospheric  argon 
cascade  arc.   These  measurements  are  made  for  a  wide 
range  of  laminar  flow  operating  conditions.   Additional 
measurements  are  made  to  determine  the  ratio  of  radiative 
to  total  wall  heat  flux  for  the  asymptotic  region  of  an 
argon  arc**.   An  important  secondary  objective  includes 


*  A  review  of  the  subject  of  anode  heat  transfer  for 
various  electrode  geometries  is  given  by  Eckert  and 

Pfender t30] . 
**  A  review  of  the  literature  pertinent  to  this  subject  is 
presented  in  the  Appendix. 
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comparison  of  the  asymptotic  wall  heat  transfer  and  electric 
field  intensity  measurements  with  the  experimental  data  of 
other  investigators  and  with  the  analytical  and  equili- 
brium and  nonequilibrium  numerical  solutions  described 
previously.   With  respect  to  the  latter,  particular  atten- 
tion is  given  to  ascertaining  both  the  validity  of  existing 
Nusselt  number  correlations  and  the  adequacy  of  single- 
fluid  models  for  the  description  of  arc  energetics.   Pur- 
suant to  gaining  additional  information  on  anode  physics 
and  the  suitability  of  anode  heat  transfer  correlations/ 
further  comparisons  are  made  with  existing  anode  heat 
transfer  models.   The  measurements  are  also  of  value  in 
providing  an  improved  understanding  of  the  development 
length  characteristics  of  the  constricted  tube,  laminar 
arc.   In  order  to  ascertain  the  influence  of  operating 
gas,  additional  experimental  results  are  presented  for 
nitrogen  and  helium  arcs. 

This  study,  which  is  the  first  serious  attempt  to 
experimentally  verify  existing  equilibrium  and  nonequili- 
brium theoretical  predictions  of  thermal  parameters  for 
the  arc  heating  region,  is  intended  to  provide  information 
useful  to  the  arc  designer  and  to  contribute  to  an  im- 
proved understanding  of  arc  physics  and  flow  modeling. 
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CHAPTER  II   ARC  GENERATOR  AND  DIAGNOSTIC  TECHNIQUES 

II. 1  Arc  Generator  and  Supporting  Systems 

A  well-defined,  stable  plasma  flow  is  achieved 

through  use  of  a  constricted  tube  plasma  arc  generator 

[341 
manufactured  by  Creare  Inc.      (Figure  II. 1).   The 

constrictor  is  formed  by  stacking  together  a  series  of 
individually  water-cooled  copper  segments  (discs)  /  each 
of  which  is  0.632  cm  wide  with  a  1.0  cm  central  bore.   Each 
segment  is  electrically  insulated  from  adjoining  segments 
by  cool  O-ring  seals  which  further  provide  for  alignment 
and  hermetic  sealing  of  the  duct.   The  lengths  of  the 
arc  heating  and  field-free  regions  may  be  varied  according 
to  the  anode  location  and  the  number  of  segments  used  to 
form  the  stack  (up  to  a  maximum  of  40) .   The  thoriated 
tungsten  cathode  is  water  cooled  and  may  be  precisely 
centered  with  respect  to  the  initial  nozzle-shaped  seg- 
ment through  micrometer  adjustment.   The  cathode-nozzle 
geometry  provides  for  a  high  degree  of  flow  symmetry  which, 
along  with  water  cooling,  provides  for  extended  cathode 
life.   If  the  tungsten  tip  of  the  cathode  is  periodically 
(every  20  hours)  re-shaped  on  a  lathe  such  that  the  proper 
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cathode-nozzle  geometry  is  maintained,  a  single  cathode 
will  provide  approximately  100  hours  of  operation. 

The  supporting  systems  which  are  necessary  to  sus- 
tain the  arc  generator  are  illustrated  schematically  in 
Figure  II. 2.   Electrical  power  is  supplied  by  a  D-C 
motor-generator  set  capable  of  delivering  75  kw  at  3  00 
amperes.   Stable  operating  conditions  are  achieved  through 
use  of  two,  step-wise  variable  ballast  resistors  which 
are  connected  in  series  with  the  power  supply  and  arc 
generator.   The  electrical  system  will  allow  for  the  use 
of  either  r-f  discharge  or  prod  starting  procedures. 

The  main  coolant  system  supplies  cooling  water  to 
the  plasma  generator  coolant  manifold  and  to  the  main 
ballast  resistor.   The  manifold  in  turn  distributes  cooling 
water  to  the  cathode,  the  starting  prod,  those  segments 
which  do  not  serve  as  calorimeters,  and  to  the  additional 
downstream  ducting  required  to  cool  the  gas  for  discharge. 
A  secondary  coolant  system  is  employed  to  reduce  flow 
demands  on  the  main  coolant  system  and  to  provide  a 
means  of  selectively  monitoring  and  controlling  the 
coolant  flow  rates  through  those  duct  segments  for  which 
wall  heat  transfer  measurements  are  desired.   This  system 
also  supplies  cooling  water  to  the  auxiliary  ballast 
resistor. 
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Receiving  its  supply  from  commerical  gas  cylinders, 
the  gas  system  facilitates  control  and  measurement  of  the 
operating  gas  flow  rate.   Stringent  purity  requirements 
are  maintained  throughout  the  gas  system  including  use  of 
research  grade  gas,  prior  ultrasonic  cleaning  of  all 
components  exposed  to  the  plasma  stream,  and  continuous 
outgassing  of  the  system  prior  to  operation. 

II. 2   Particular  Operating  Procedures 
and  Limitations 

Despite  the  aforementioned  purity  requirements,  the 
level  of  duct  contamination  resulting  from  arc  operation 
has  steadily  increased  over  the  past  two  year  period. 
During  the  earlier  stages  of  this  research,  when  the  major 
portions  of  the  electric  field  intensity  and  total  wall 
heat  flux  data  were  acquired,  there  were  no  visible  indi- 
cations of  surface  contamination,  even  after  extended  per- 
iods of  operation  (up  to  40  hours).   However,  during  the 
latter  phases  of  the  research,  when  the  radiative  wall 
heat  flux  data  were  acquired,  the  level  of  contamination 
continually  increased  in  spite  of  special  precautions  that 
were  taken  in  addition  to  those  listed  above. 

The  contamination  is  evidenced  by  a  thin  deposit 
which  coats  the  duct  walls  in  a  distinct  and  repetitive 
pattern.   The  deposit,  which  appears  to  be  a  dark  oxide 
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layer,  forms  on  the  leading  (upstream)  edge  of  each  of  the 
segments  in  the  arc  heating  region  and  diminishes  in  the 
axial  direction,  leaving  the  downstream  half  of  the  seg- 
ment surfaces  free  of  visible  contamination.  The  deposit, 
which  at  times  has  a  bluish  cast,  has  not  been  observed  in 
the  field-free  region,  not  even  on  the  segment  immediately 
downstream  of  the  anode. 

Similar  deposits  have  been  observed  by  other  in- 

T13  32  331 
vestigators    '       '         ;  however  opinions  differ  concerning 

the  effect  of  surface  contamination  on  experimental  results. 

r  13] 

Emmons     claims  that  the  overall  arc  characteristics  are 

virtually  unaffected  by  small  amounts  of  contamination. 

[33] 
On  the  other  hand,  Runstadler     reports  a  reduction  in 

the  electric  field  intensity  for  helium  of  some  10%  when 

data  are  acquired  under  contamination-free  conditions. 

Throughout  the  course  of  this  study,  both  electric  field 

intensity  and  total  wall  heat  flux  measurements  have  been 

consistently  reproducible  to  within  their  respective 

experimental  uncertainties,  irrespective  of  the  level  of 

contamination.   None  the  less,  it  was  assumed  initially 

that  contamination  would  render  the  proposed  radiation 

gage  completely  useless,  since  the  dark  deposits  would 

significantly  alter  the  optical  properties  of  the  gage 

surface.   While  useful  data  were  acquired  from  the 
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radiation  gage,  the  contamination  level  within  the  duct 
still  prevented  the  gage  from  being  used  as  it  was 

originally  intended  (see  the  Appendix) . 

r  30 1 

As  reported  by  McKee /  et  al.     ,  the  difficulties 

associated  with  obtaining  reproducible  measurements  of 
anode  characteristics  are  most  likely  due  to  the  pro- 
nounced affect  that  surface  contamination  has  on  electrode 
performance.   In  connection  with  this  problem,  a  number 
of  investigators    '        have  devoted  considerable 
effort  to  the  elimination  and  investigation  of  surface 

contamination  within  constricted  tube  arcs.   In  summarizing 

[33] 
work  performed  in  this  area  at  Dartmouth,  Runstadler 

suggests  a  number  of  procedures  which  may  be  used  to  re- 
duce the  level  of  contamination  within  the  duct.   These 
recommended  procedures  include  use  of  ultra-pure  grade 
gas,  elimination  of  possible  cathode  contaminants  through 
use  of  a  sonic  orifice,  outgassing  of  the  duct  surfaces, 
and  nickel  plating  of  all  copper  surfaces  exposed  to  the 
plasma  stream.   The  plating  procedure  was  found  necessary 
in  order  to  eliminate  outgassing  from  the  copper  surfaces 
during  arc  operation. 

Similar,  yet  considerably  less  involved,  measures 
were  taken  in  this  study  in  an  attempt  to  eliminate  the 
contamination  of  the  duct  surfaces.   Ultra-pure  grade 
(as  opposed  to  commercial  grade)  gas  was  used;  however 
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its  use  did  not  result  in  a  significant  reduction  in 
surface  contamination.   This  is  not  to  say,  however,  that 
use  of  commercial  grade  gas  will  not  result  in  surface 
contamination.   After  a  thorough  cleaning  of  the  gas 
system  failed  to  effect  any  change  in  the  degree  of  sur- 
face contamination,  the  inlet  portion  of  the  gas  system 
was  virtually  eliminated  by  connecting  the  gas  cylinder 
directly  to  the  plasma  generator  (a  short  piece  of  plastic 
tubing  was  connected  between  the  gas  inlet  to  the  plasma 
generator  and  the  gas  regulator  and  the  regulator  was 
connected  directly  to  the  gas  cylinder).   In  addition, 
the  gas  regulator  (a  two-stage,  metal  diaphram  type)  was 
completely  disassembled,  cleaned  ultra-sonically  with 

acetone  and  isopropyl  alcohol,  and  allowed  to  degas  at 

—  6 
10   mm  Hg  for  a  period  of  24  hrs.   Several  hours  prior 

to  starting  the  plasma  generator,  the  combined  generator 
and  gas  system  was  evacuated  to  approximately  40  microns 
and  the  metal  components,  particularly  the  duct  segments, 
were  heated  to  accelerate  degassing.   Unfortunately,  none 
of,  or  any  combination  of,  these  procedures  led  to  any 
significant  reduction  in  the  level  of  surface  contamina- 
tion.  In  addition,  the  possible  influence  of  prod  and 
r-f  discharge  starting  techniques  was  investigated,  as 
were  the  procedures  used  to  mechanically  and  chemically 
clean  the  duct  segments.   Special  care  was  taken  in 
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positioning  the  cathode  within  the  initial  nozzle  segment 
so  as  to  minimize  possible  contamination  from  erosion  of 
the  tungsten  cathode  tip.   Again,  there  were  no  per- 
ceptable  changes  in  the  level  of  surface  contamination  with 
any  of  these  additional  measures. 

On  the  basis  of  these  results,  it  is  tentatively 
concluded  that  the  surface  contamination  is  an  oxide 
layer  which  is  formed  by  a  chemical  reaction  between  the 
copper  duct  walls  and  the  gases  which  escape  from  these 

walls  during  operation  of  the  arc.   This  conclusion  is 

T331 
supported  by  the  observations  of  Runstadler     concerning 

the  nickel  plating  of  the  copper  surfaces,  and  also  by 

the  experiments  performed  in  this  study  with  the  thin 

film  aluminum  surface  on  the  radiation  gage  (Appendix) . 

Further  support  of  this  conclusion  is  provided  by  the 

[321 
x-ray  diffraction  analysis  conducted  by  McKee ,  et  al. 

which  identified  the  surface  contaminants  as  copper-oxide 
and  carbon.   While  not  fully  understood,  the  form  of  the 
axial  distribution  of  the  surface  contamination  on  seg- 
ments in  the  arc  heating  region  and  the  absence  of  surface 
contamination  in  the  field-free  region  indicate   that 
certain  electrical  effects  must  also  influence  the 
contamination . 
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In  addition  to  experiencing  contamination  problems, 
some  difficulties  arose  in  connection  with  generating  He 
and  N_  arc  plasmas.   The  maximum  current-carrying  capa- 
bilities of  the  starting  prod  (75  amp)  and  r-f  discharge 
(50  amp)  circuits,  while  adequate  for  initiating  argon 
arcs,  are  insufficient  to  permit  starting  in  either 
nitrogen  or  helium.   Accordingly,  a  special  pre-start 
procedure  was  devised.   This  procedure  involves  use  of  a 
dual  gas  supply  manifold  whereby  both  Ar  and  either  N2  or 
He  can  be  supplied  to  the  arc  generator  simultaneously. 
Following  initiation  of  a  stable,  50  amp  Ar  arc,  the 
proportion  of  N?  or  He  to  Ar  is  slowly  increased  from  0  to 
100%  while  the  open-circuit  voltage  is  steadily  increased 
to  maintain  a  stable  operating  condition.   While  this 
procedure  is  adequate  for  initiating  N2  and  He  arcs ,  the 
maximum  obtainable  open-circuit  voltage  of  the  power 
supply  is  nonetheless  insufficient  to  sustain  stable 
operation  of  these  gases   over  as  wide  a  range  as  is 
possible  with  Ar. 

The  restrictions  which  are  placed  upon  both  high 
and  low  current  operation  of  the  arc  generator  are  a  re- 
sult of  certain  physical  limitations  of  the  power  supply 
and  external  ballast  resistor.   Specifically,  high  current 
operation  is  limited  by  the  maximum  open-circuit  voltage 
capability  of  the  power  supply.   In  contrast,  low  current 
operation  is  restricted  by  the  lower  limit  of  electrical 
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stability  which,  in  turn,  is  a  function  of  the  available 
ballast  resistance  as  well  as  the  maximum  obtainable  open- 
circuit  voltage  of  the  power  supply. 

These  limits  of  operation  can  be  explained  by  con- 
sidering a  set  of  hypothetical  (yet  typical)  load  line 
characteristics  for  an  arc  generator  and  its  associated 
D-C  power  supply  (Figure  II. 3).   As  indicated  by  the 
electrical  circuit  which  accompanies  Figure  II. 3,  the 
power  supply  load  line  characteristic  is  a  plot  of  the 
voltage  available  to  the  arc  generator  versus  the  available 
arc  current.   Generally,  the  power  supply  characteristic  is 
represented  by  a  straight  line  which  intersects  the  vol- 
tage axis  at  the  value  of  the  open-circuit  voltage  and 
which  has  a  slope  equal  to  the  negative  of  the  sum  of 
the  constant  internal  resistance  of  the  power  supply  (R  ) 

S3 

and  the  variable  external  ballast  resistance  (R,  )  .   The 

b 

arc  generator  load  line,  on  the  other  hand,  represents  a 
plot  of  the  current  conducted  by  the  arc  versus  the 
voltage  required  to  sustain  that  current.   On  the  whole, 
the  arc  generator  is  characterized  by  a  wide  spectrum  of 
load  lines  which  result  from  various  combinations  of  the 
operating  gas,  the  operating  gas  flow  rate  (m) ,  the 
operating  pressure  (p) ,  the  overall  (cathode  to  anode)  duct 
length  (L) ,  and  the  duct  radius  (R) . 
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The  actual  operating  state  (for  a  specified  gas, 
A,  p,  L,  and  I)  of  a  constant  diameter  arc  generator  is 
uniquely  determined  by  the  point  of  intersection  of  the 
load  (arc  generator)  and  supply  characteristics.  The 
degree  of  electrical  stability  which  exists  for  a  given 
operating  state  is  a  function  of  the  relative  slopes  of 
the  two  load  lines  at  their  point  of  intersection.  The 
criterion  that  must  be  satisfied  in  order  to  achieve 
stable  operation  is  given  by 


<fl>         -  <SS>  -  o  (ii. i) 

LOAD  SUPPLY 


which,  in  the  case  of  a  linear  supply  characteristic, 
reduces  to 


(S)      +  R+  R,  -  0  (II. 2) 

51  LOAD     S     b 


It  is  noted  at  this  point  that,  while  Equation  II. 2 
establishes,  in  principle,  the  criterion  for  stable  opera- 
tion, in  practice  the  slope  of  the  load  characteristic  must 
always  be  greater  than  the  slope  of  the  supply  character- 
istic if  stable  operation  is  to  be  realized.   That  is, 
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dF 
(~)      +  R   +  R,   >>  0.  (II. 3) 

dI  LOAD     S     b 


dE 
Since  R   is  constant  and  since  {t=-)  is  dictated  by  the 

LOAD 
desired  arc  operating  condition,  improved  electrical  sta- 
bility at  each  condition  can  be  effected  only  by  increasing 
R,  .   This,  in  turn,  necessitates  an  increase  in  the  open- 
circuit  voltage.   Accordingly,  for  stable  arc  operation, 
it  is  generally  advisable  to  use  the  highest  possible 
settings  of  both  the  ballast  resistor  and  the  open-circuit 
voltage  of  the  power  supply. 

Irrespective  of  the  value  of  R,  ,  Equation  II. 3  will 

be  satisfied  as  long  as  operation  is  confined  to  the 

* 
positive  slope  portions  of  the  arc  generator  load  lines  . 

Nevertheless,  the  maximum  open-circuit  voltage  capability 
of  the  power  supply  can  place  a  restriction  on  high  cur- 
rent operations.   For  example,  if  all  of  the  ballast 
resistance  has  been  removed  in  order  to  operate  in  condi- 
tions (Figure  II.  3),  further  operation  at  higher  arc 
currents  (position-a ' )  or  operation  at  greater  values  of 
m  and/or  L  (position-a")  can  be  achieved  only  by  increasing 


*  Note  that  operation  at  the  initial  point  of  intersection 
(e.g.,  points  x  and  y  in  Figure  II. 3)  is  not  likely  since, 
according  to  Equation  II. 2,  the  arc  would  be  inherently 
unstable  at  this  condition. 
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the  open-circuit  voltage.   Therefore,  aside  from  possible 
over  heating  of  the  electrodes,  high  current  operation 
is  limited  by  the  maximum  obtainable  open-circuit  voltage 
of  the  power  supply. 

According  to  the  stability  requirement  expressed  by 
Equation  II. 3,  operation  on  the  negative  slope  portions  of 
the  arc  generator  load  lines  (e.g.,  position-b)  is  in- 
herently less  stable  than  operation  in  the  positive  slope 
regions.   In  addition,  a  lower  limit  for  stable  operation 
exists  at  the  current  for  which  the  slopes  of  the  load  and 
supply  characteristics  are  equal  (e.g.,  position-c) .   In 
practice,  the  arc  will  extinguish  prior  to  this  lower 
limit;  however,  lower  current  operation  (position-c1)  or 
operation  at  higher  valves  of  m  and/or  L  (position-c")  can 
be  achieved  by  increasing  the  open-circuit  voltage  and 
increasing  R,  .   This  procedure  can  be  continued  until 
the  maximum  open-circuit  voltage  of  the  power  supply  (or 
the  maximum  value  of  R,  )  is  reached. 

For  the  duct  lengths  and  operating  gas  flow  rates 


employed  in  this  study,  the  point  at  which  the  slope  of 

dE. 
dl 


dE 
the  load  characteristic  is  equal  to  zero  (  (^f)      =  0) 

LOAD 


occurs  at  an  arc  current  of  approximately  25  amp  for  Ar , 
150  amp  for  N  ,  and  at  a  current  in  excess  of  400  amp  in 
He.  Therefore,  according  to  the  previous  arguments,  low 
current  operations  in  all  three  gases  and  high  current 
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operations  in  N2  and  He  are  limited  by  the  power  supply. 
High  current  operation  in  Ar  is  limited  by  the  maximum 
power  that  can  be  dissipated  by  the  main  ballast  resistor. 
It  is  further  noted  that  over  heating  of  the  anode  (as 
evidenced  by  ablation  of  the  anode  surface)  occurred  at 
the  higher  arc  currents  in  N~  (I  >  200  amp)  and  in  He 
(I  >  100  amp) . 

II. 3  Diagnostic  Methods  and  Related  Instrumentation 
The  overall  and  local  arc-column  characteristics 
are  determined  by  measuring  the  potential  difference 
between  the  "floating"  duct  segments  and  the  cathode.   The 
arc  voltage  measurements  are  used  both  to  determine  the 
axial  distribution  of  the  electric  field  intensity  and, 
along  with  the  arc  current  and  calorimetric  measurements, 
to  determine  the  mixed-mean  enthalpy  of  the  gas  at  each 
segment  location.   The  arc  voltages  are  measured  with  a 
high  input  impedance  (4  megaohm)  strip-chart  recorder 
(Electronik-19  by  Honeywell) .   The  arc  current  is  deter- 
mined by  monitoring  the  voltage  drop  across  a  calibrated 
shunt  resistor  with  a  digital  voltmeter  (Dana,  model 
#  54030).   As  shown  schematically  in  Figure  II.  4,  further 
circuitry  is  available  to  permit  up  to  400  ma  of  current 
to  be  drawn  from  a  duct  segment  while  its  potential  (with 
respect  to  the  cathode)  is  monitored. 
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Figure  II. 4   Circuit  Used  for  Impressed-Current  (Ia) 
Method  of  Wall  Potential  Measurement 
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The  energy  loss  to  the  duct  segments/  to  the 
cathode,  and  to  the  remaining  arc  generator  components  is 
determined  calorimetrically .   Accordingly,  each  of  these 
components  is  equipped  with  a  specially  designed  tempera- 
ture difference  (AT)  transducer.   The  AT  transducers 
(Figure  II. 5)  employ  a  single  pair  of  copper-constantan 
thermocouples,  mounted  in  plexiglass  inlet  and  outlet 
plenum  chambers,  and  wired  in  series  to  provide  a  voltage 
signal  which  is  directly  proportional  to  the  coolant 
temperature  rise.   Considerable  care  was  taken  in  the 
shielding  and  calibration  (against  mercury-in-glass 
thermometers)  of  these  transducers  and  they  are  accurate 
to  within  5%  of  the  indicated  reading.   The  transducer 
output  voltages  are  amplified  (Dana  model  #3850-V2)  prior 
to  being  displayed  on  the  digital  voltmeter.   The  abso- 
lute temperature  of  the  operating  gas  at  both  inlet  and 
discharge  locations  is  measured  with  a  copper-constantan 
thermocouple  (as  referenced  to  a  copper-constantan  ice- 
junction)  .   The  gas  temperatures  are  monitored  by  the 
digital  voltmeter. 

A  unique  radiative  heat  transfer  gage  is  used  to 
delineate  the  radiative  contribution  to  the  total  wall 
heat  flux.   As  shown  in  Figure  II. 6,  the  radiation  gage 
consists  of  two  semicylindrical  copper  segments  which 
are  thermally  isolated  from  one  another,  individually 
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water  cooled,  and  supplied  with  individual  AT-transducers. 
The  inside  wall  of  each  half  segment  is  coated  with  either 
a  highly  reflecting  (aluminum)  or  highly  absorbing  (black) 
surface.   The  overall  dimensions  of  the  assembled  gage 
are  precisely  the  same  as  those  of  a  standard  duct  seg- 
ment.  A  thorough  discussion  of  the  principle  of  opera- 
tion, the  fabrication,  and  the  associated  experimental 
problems  involved  in  the  use  of  the  radiation  gage  are 
left  to  the  Appendix. 

The  operating  gas  flow  rate  and  the  coolant  flow- 
rates  through  the  cathode  and  those  duct  segments  serving 
as  calorimeters  are  measured  with  precision  accuracy 
(±  1%)  flowmeters  (Fischer  &  Porter  model  no.  10A0735) . 
The  two  gas  flowmeters  have  a  combined  range  of  from  0.02 
to  1.50  lb/min  in  argon.   At  the  time  of  each  flowmeter 
reading,  the  inlet  gas  temperature  and  pressure  are 
measured  and  used  to  correct  the  indicated  flow  rate 
should  they  differ  from  the  temperature  and  pressure  of 
calibration.   An  additional  precision  accuracy  flowmeter 
is  used  to  measure  the  total  coolant  flow  rate  to  the 
plasma  generator. 

With  measurement  of  the  arc  current  (I)  and  axial 
voltage  distribution  (E) ,  the  gas  inlet  temperature  (TINLET) 
and  flow  rate  (m) ,  and  the  power  losses  to  the  duct 
segments  (PLS)  and  to  the  cathode  (PLC) ,  the  axial 
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distribution  of  the  gas  enthalpy  can  be  determined  by 
performing  energy  balances  at  successive  axial  locations. 


) 
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Figure  II. 7  Energy  Balance  Diagram 


With  reference  to  Figure  II. 7,  the  enthalpy  of  the  gas 
at  the  exit  plane  of  a  segment  in  the  arc  heating  region 
is  given  by 


hj  =  Cp  TINLET  + 


E 


E.I-  Z    (PLS)±  -  PLC 


(II. 3) 


m 
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This  value  is  then  algebraically  averaged  with  the  enthalpy 
determined  by  the  previous  (j-1)  energy  balance,  and 
the  resulting  "mixed-mean"  enthalpy  is  assigned  to  an 
axial  location  corresponding  to  the  center  of  the  j 
segment. 

To  experimentally  verify  the  accuracy  of  the  diag- 
nostic procedures  employed  in  this  study,  an  attempt  was 
made  to  "close  off"  the  energy  balance  for  the  entire 
plasma  generator.   While  the  precision  of  this  energy 
balance  is  limited  by  unavoidable  ambient  losses,  it  was 
still  possible  to  obtain  closure  to  within  the  cumulative 
experimental  uncertainty  of  the  component  measurements 
(±  20%) .   A  more  precise  check  is  one  based  on  a  local 

energy  balance  performed  on  a  segment  located  in  the 

* 
fully-developed  portion  of  the  arc  heating  region  .   In 

every  case,  closure  of  this  local  energy  balance  was 
obtained  to  within  5%  -  a  value  comparable  to  the  total 
uncertainty  of  these  measurements  (±5-10%) .   The  consistent 
reproducibility  of  data  taken  over  a  two  year  period  pro- 
vides further  confidence  in  the  accuracy  of  the  measure- 
ments as  well  as  the  consistency  of  the  overall  system. 

A  detailed  uncertainty  analysis  is  used  for  estima- 
ting the  uncertainties  associated  with  each  measured 


*  The  nature  of  this  energy  balance  will  be  discussed 
in  detail  in  the  following  chapter  (Section  III.  3). 
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variable  and  for  calculating  the  propagation  of  these 
uncertainties  into  each  reported  result.   The  method  used 
in  this  analysis  is  the  standard  approach  for  single-sample 

r  36 1 

experiments     .   The  uncertainty  calculations,  as  well 
as  the  conversion  and  reduction  of  raw  data,  are  per- 
formed on  the  computer.   The  thermo-physical  properties 

required  for  determining  the  mean  temperature.  T  =T  (h  ) , 

^  c  m  m  m 

and  the  Nusselt  number-Ohmic  heating  parameter  correlations 

[371 
(k  and  a  )  are  obtained  from  DeVoto     and  Drellishak, 
m      m 

T  3  8 1  T  3  9  ] 

Knopp,  and  Cambel     in  the  case  of  Ar ;  Lick  and  Emmons 

and  DeVoto  and  Li     in  the  case  of  He;  and  Drellishak, 

[41]  r  4  2 1 

Aeschliman,  and  Cambel  and  Avco  for  N2»  The  re- 
ported uncertainties  of  these  properties  are  also  consi- 
dered in  the  aforementioned  uncertainty  analysis.   A  more 

detailed  description  of  the  arc  facility,  diagnostic  tech- 

f  43] 
niques,  and  uncertainty  analysis  is  given  by  Lukens 


37 


CHAPTER  III   EXPERIMENTAL  RESULTS 

111,1   Range  of  Measurements 
Experimental  data  were  acquired  for  Ar ,  N~,  and 
He  arcs  over  the  following  ranges  of  laminar  flow  operating 
conditions : 

Argon;     I  =  35  -  2  40  amp,  m  =  0.03  -  0.11  lb/min , 
constrictor  length  =  5  in  (20  segments) . 

Nitrogen;  I  =  100  -  200  amp,  m  =  0.03  lb/min , 

constrictor  length  =  5  in  (20  segments) . 

Helium;    I  =  50  amp,  rfi  =  0.07  lb/min, 

constrictor  length  =  2.5  in  (10  segments). 

All  measurements  were  made  with  a  1.0  diameter  copper 
duct  operating  at  near  atmospheric  pressure.   The  assertion 
of  laminar  flow  is  based  upon  the  flow  transition  studies 
by  Incropera  and  Leppert     and  by  Runstadler     in 
which  the  onset  of  turbulence  is  reported  to  occur  at  a 
mean  Reynolds  number  (based  on  properties  evaluated  at  T  ) 
of  700  to  900.   The  largest  value  of  the  mean  Reynolds 
number  encountered  in  this  study  is  approximately  600 
(Argon,  I  =  35  amp,  m  =  0.11  lb/min).   Hence,  the  above 
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stated  ranges  of  operating  conditions  are  assumed  to  be 
representative  of  laminar  flow. 

Rather  than  presenting  a  complete  tabulation  of  the 
data  acquired  for  the  inlet  region,  certain  data  were 
selected  to  demonstrate  the  effect  of  varying  gas  flow 
rate  and  arc  current  on  the  development  characteristics 
of  the  arc.   All  of  the  acquired  data  pertaining  to  the 
asymptotic  region  are  presented  along  with  related  data 
from  other  investigators.   The  interpretation  of  the 
experimental  results  as  well  as  comparisons  with  related 
theories  will  be  the  subject  of  Chapter  IV. 

III. 2   Electrical  Characteristics 
Prior  to  presenting  the  experimental  data  for  the 
arc-column  (electrical)  characteristics/  a  discussion  is 
in  order  concerning  the  method  employed  in  the  acquisition 
of  this  data. 

A  method  frequently  used  in  the  measurement  of 
arc-column  characteristics  involves  use  of  the  duct  seg- 
ments as  "floating"  potential  probes [1°  '12  '15  ,32] .   It 
is  commonly  recognized,  however,  that  a  floating  duct  wall 
will  assume  a  potential  which  is  negative  with  respect 

to  the  plasma  core  by  an  amount  equivalent  to  the  plasma 

* 
sheath  potential  .   Accordingly,  on  the  assumption  that 


*  Further  treatment  on  the  formulation  and  characteristics 
of  plasma  sheaths  is  given  in  references  [30]  and  [45] . 
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the  sheath  potential  is  appreciable,  some  investigators 
have  questioned  the  suitability  of  using  a  floating  duct 
wall  to  measure  plasma  potentials. 

In  reference  to  their  cascade  arc  potential  measure- 
ments, this  question  has  been  raised  most  strongly  by 

r  6  71 

Pf ender ,  et  al .      .   In  performing  their  measurements, 
they  allowed  a  current  (I  )  to  be  drawn  from  the  duct 
wall  while  they  simultaneously  monitored  the  potential  of 
this  wall  with  respect  to  the  cathode.   These  measure- 
ments indicated  that  the  recorded  wall  potential  was 
strongly  influenced  by  the  magnitude  of  the  current  drawn. 
In  particular,  they  observed  that  the  voltage  increases 

by  40  to  100%  over  the  floating  wall  (I   =  0)  potential 

a 

when  a  current  of  up  to  50  ma  is  drawn.   Above  50  ma,  the 
segment  potentials  were  observed  to  be  independent  of  the 
current  drawn.   With  regard  to  this  observation,  they 
argued  that,  for  I   >  50  ma,  the  charge  separation  which 
otherwise  would  exists  in  the  plasma  sheath  region  is 
exactly  cancelled  and,  consequently,  the  duct  wall  is  at 
the  potential  of  the  plasma  column.   They  therefore  con- 
cluded that  floating  wall  potential  measurements  do  not 
provide  a  suitable  indication  of  the  actual  plasma 
potential,  but  that  accurate  measurements  of  this  para- 
meter can  be  made  by  recording  the  wall  potential  while 
a  current  in  excess  of  50  ma  is  drawn  from  the  wall. 
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[12  131 
In  contrast/  Emmons    '     indicates  that  the  sheath 

potentials  for  an  argon-copper  duct  combination  are 

approximately  1  to  2  volts,  for  which  case  the  floating 

wall  potentials  are  suitably  close  to  those  in  the  arc 

column.   In  support  of  this  contention  is  the  experimental 

technique  used  by  Maecker     wherein  every  attempt  was  made 

to  eliminate  impressed  wall  currents  (I  )  while  performing 

a 

*  T461 

plasma  potential  measurements  .   Schreiber     ,  on  the 

other  hand,  has  measured  significant  voltage  increases 
with  current  drawn,  but  he  attributes  these  increases  to 
voltage  drops  across  the  relatively  cool  gas  boundary  lay- 
ers.  Additional  support  for  Schreiber' s  deduction  can 

be  inferred  from  the  measurements  made  by  researchers  at 

[32  331 
Dartmouth    '    .   While  the  Dartmouth  data  were  used 

primarily  in  support  of  anode  studies,  the  arc  configura- 
tion and  external  circuitry  used  for  these  measurements 

are  similar  to  those  employed  in  this  study  (Figure  II. 4) 

[71 

and  to  those  used  by  Pfender,  et  al.    . 

In  an  attempt  to  resolve  these  inconsistencies,  the 
variation  in  the  duct  wall  potential  with  current  drawn 
was  measured  at  various  axial  locations  and  arc  currents. 
As  shown  in  Figure  III.l,  these  measurements  indicate  that 


*  Maecker  contends  that  the  measured  wall  potential 

difference  between  two  segments  will  be  in  error  by  an 

amount  proportional  to  the  impressed  current,  I  . 

a 
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the  percentage  increase  in  wall  potential  with  increasing 

I  diminishes  both  with  increasing  axial  distance  and 
a 

with  increasing  arc  current.   They  also  indicate  that 

the  wall  potential  increases  continuously  with  increasing 

I   over  the  entire  range  of  I   values  investigated.   In 
a  a 

this  respect,  the  distributions  of  Figure  III.l  are 


qualitatively  similar  to  the  potential  distributions 

urre 
[33] 


[151 
recorded  by  Maecker     and  to  the  lower  current  ranges 


of  anode  conduction  reported  by  Runstadler 

Simultaneous  calorimetric  and  wall  potential  mea- 
surements indicate  that  the  total  wall  heat  flux  to  a 
duct  segment  does  not  vary  appreciably  (<  3%)  over  the 

entire  range  of  I  .   In  addition,  the  value  of  the 

a 

asymptotic  electric  field  intensity  (E    )  was  observed 

x ,  a 

to  be  independent  of  the  current  drawn.   This  fact 

is  discussed  by  Emmons  and  Land     and  by  McKee ,  et  al. 

and  is  strongly  supported  by  the  experimental  measure- 

[141 
ments  of  Morris,  et  al.     ,      wherein  the  voltage  gradient 

was  measured  using  three  independent  methods.   With  the 

first  method,  the  value  of  E    was  computed  from  the 

x ,  a 

asymptotic  energy  balance  equation  (this  equation  is  dis- 
cussed later  in  this  chapter)  along  with  calorimetric 
measurements  of  the  power  losses  to  the  duct  segments. 
The  second  method,  which  is  similar  to  that  employed  by 
Maecker,  utilizes  two  of  the  duct  segments  as  potential 
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probes.   Unlike  Maecker ' s  method,  however,  a  small  current 
was  supplied  to  the  arc  through  these  probes  (in  an  attempt 
to  cancel  the  voltage  drop  across  the  "contact  resistance" 
between  the  probes  and  the  arc  column)  while  the  potential 
difference  between  the  two  duct  segments  was  measured. 
This  potential  drop  was  then  divided  by  the  distance  be- 
tween the  two  probe  segments  to  give  the  value  of  E    .   The 

x ,  a 

third  method  used  by  Morris,  et  al .  involved  a  cathode 
withdrawal  experiment  wherein  the  total  arc  current  was 
held  constant  while  the  total  (cathode  to  anode)  voltage 
drop  was  recorded  as  a  function  of  the  total  (cathode 
to  anode)  duct  length.   These  three  methods,  when  applied 
to  the  central  region  of  a  hydrogen  arc,  resulted  in  values 
of  the  electric  field  intensity  which  agreed  to  within 
approximately  10%  -  a  value  comparable  to  the  combined 
uncertainties  of  these  measurements. 

The  agreement  demonstrated  by  these  experiments  is 
presumably  a  result  of  the  conditions  which  characterize 
the  fully-developed  region.   That  is,  regardless  of  the 
magnitude  of  possible  radial  electric  field  components, 
the  voltage  difference  between  any  two  axial  points  in 
the  fully-developed  region  is  invariant  as  long  as  the  two 
points  are  at  the  same  radial  location.   If  this  were  not 
the  case,  the  axial  component  of  the  electric  field 
intensity  would  vary  with  axial  location  and  the  flow 
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would  no  longer  be  asymptotic.   The  independence  of  E 

with  I   was  tested  under  the  most  severe  condition 
a 

(I  =  50  amps,  I   =  400  ma)  of  this  investigation  and,  as 
a 

shown  in  Figure  III. 2,  the  slopes  of  the  linear  portions 

of  the  wall  potential  distributions  (E    )  are  the  same 

x ,  a 

for  I   =0  and  I   =  400  ma.   It  is  further  noted  that  the 
a  a 

electric  field  intensity  distribution  on  either  side  of  a 

segment  from  which  current  is  being  drawn  is  unaffected 

by  the  wall  potential  variations  of  that  segment. 

Although  it  is  concluded  that  the  current  drawn 

from  a  duct  segment  (up  to  I   =  400  ma)  does  not  apprecia- 

a 

bly  effect  the  total  wall  heat  flux  or  E    ,  the  increase 
1  x,a 

in  wall  potential  with  increasing  I   will  nonetheless  have 
a  direct  effect  on  the  energy  balance  measurements  used 
to  determine  the  mean  enthalpy  of  the  gas,  particularly 
at  the  lower  arc  currents.   As  indicated  in  Figure  III. 3, 
the  mean  enthalpy  of  a  50  amp,  argon  plasma  is  considerably 
greater  when  computed  using  wall  potentials  measured  with 
I   =  400  ma  than  it  is  when  computed  with  floating  wall 
potentials.   While  the  relative  accuracy  of  these  two 
distributions  could  best  be  determined  by  comparison  with 
experimental  enthalpy  profile  data,  such  data  are  not 

available.   However,  a  comparison  is  made  with  the  non- 

r  25 1 

equilibrium  numerical  solution  by  Clark     .   This 

solution,  as  applied  to  the  asymptotic  region  of  an 
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atmospheric  argon  plasma  in  a  1.0  cm  diameter  duct,  pre- 
dicts electric  field  intensity  and  total  wall  heat  flux 

values  which  agree  with  measured  values  to  within  the 

* 

experimental  uncertainties  (<  6%)  .   As  indicated  in 

Figure  III. 3,  this  solution  also  predicts  an  asymptotic 
mean  enthalpy  which  agrees  well  with  the  experimental 
value  determined  from  the  use  of  floating  wall  potentials 
While  the  conclusions  that  can  be  drawn  from  this  compari- 
son are,  at  best,  tenuous,  the  agreement  between  theory 
and  experiment  shown  in  Figure  III.  3  indicates  that  the 
floating  duct  wall  serves  as  the  more  suitable  potential 
probe . 

After  modifying  the  wall  boundary  conditions  to 

include  current  drawn  from  the  duct  wall,  the  nonequili- 

T251 
brium  solution  by  Clark      (which  includes  a  plasma 

sheath  model)  was  used  to  investigate  the  effect  of  I 

a 

on  the  plasma  sheath  potential.   At  I  =  50  amp  and 

I  =  400  ma,  the  nonequilibrium  solution  predicts  a 
a 

decrease  (from  its  floating  wall  value)  in  the  plasma 
sheath  potential  from  3.60  volts  to  1.75  volts.   This 
predicted  reduction  in  sheath  potential  is  considerably 
less  than  the  increases  in  wall  potential  measured  under 


*  Comparisons  between  experiment  and  the  nonequilibrium 
model  are  presented  in  Chapter  IV. 
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* 

the  same  conditions  . 

On  the  basis  of  the  measured  wall  potential  dis- 
tributions of  Figure  III.l,  and  on  the  basis  of  the 
comparisons  with  the  nonequilibrium  solution,  it  is  con- 
cluded that  the  observed  increases  in  wall  potential  are 
primarily  a  result  of  the  wall  current  (I  )  being  con- 

ducted  across  the  thin  layers  of  relatively  high  electri- 
ze* 
cal  resistance  gas  adjacent  to  the  duct  walls   .   Accord- 
ingly, it  is  felt  that  the  "floating"  wall  potential 
provides  a  more  accurate  measure  of  the  actual  plasma 
column  potential,  at  least  for  those  conditions  where  the 
plasma  sheath  potential  is  small  in  comparison  to  the 
wall  potential  (on  the  basis  of  the  nonequilibrium  solu- 
tion, the  sheath  potential,  which  increases  with  in- 
creasing arc  current,  attains  a  maximum  value  of  approx- 
imately 5.0  volts  at  I  =  240  amp).   Therefore,  all  of 
the  following  data  pertaining  to  the  electrical  character- 
istics of  the  arc-column  are  based  upon  floating  wall 
potential  measurements. 

The  axial  distribution  of  the  wall  potential  for 
an  argon  arc  and  its  variation  with  arc  current  and  with 


*  A  decrease  in  the  plasma  sheath  potential  results  in  an 
increase  in  the  duct  wall  potential  measured  with 
respect  to  the  cathode. 

**This  process  is  treated  in  greater  detail  by 
Runstadler ^333 . 
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operating  gas  flow  rate  are  shown  in  Figures  III.  4  and 
III. 5,  respectively.   Aside  from  the  anode  fall  region,  the 
wall  potential  varies  linearly  with  axial  distance  for 
entry  lengths  greater  than  2.5  inches.   The  0.25  in 
(1  segment  width)  location  of  the  zero  potential  intercept 
is  due  to  the  positioning  of  the  cathode  within  the  initial 
nozzle-shaped  segment. 

As  indicated  in  Figure  III.  5,  the  increase  in 
potential  with  increasing  flow  rate  is  to  be  expected 
since,  with  an  increase  in  the  thermal  development  length 
(and  thus  an  increase  in  the  overall  arc  resistance)  ,  a 
larger  potential  must  be  applied  to  maintain  a  constant 
arc  current.   In  contrast,  if  the  fully-developed 
temperature  profiles  are,  in  fact,  independent  of  gas  flow 
rate  (as  they  must  be  in  the  presence  of  negligible  viscous 
dissipation) ,  then  the  asymptotic  value  of  the  electric 
field  intensity  (E    )  should  not  vary  with  gas  flow 

X  a  EL 

rate.   Such  is  the  case  for  the  lower  flow  rates  shown 

in  Figure  III.  5.   The  slight  variation  in  E    with  higher 

x ,  a 

flow  rates  is  attributed  to  possible  departures  from 
axisymmetric  flow  and  nonuniform  arc  attachment  at  the 
anode  surface. 

The  asymptotic  values  of  the  electric  field  inten- 
sity for  Ar ,  N  ,  and  He  and  their  dependence  on  arc 
current  are  plotted  in  Figure  III. 6.   Included  for  compari- 
son are  the  experimental  distributions  of  Runstadler 
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and  Emmons  and  Land      in  the  case  of  Ar ,  Emmons 


[321  . 
and  McKee ,  et  al .      in  the  case  of  He,  and  Emmons 


[13] 


[13] 
for  N2.   While  Emmons     contends  that  there  are  no 

regions  (in  a  tube  with  L/D  =  14)  for  which  his  N?  data 
are  independent  of  laminar  flow  rate,  asymptotic  condi- 
tions (as  evidenced  by  constant  values  of  E   and  q  in  the 
axial  direction)  were  observed  in  this  study  for  a 
L/D  -   12.6.   Consequently,  Emmons'  lowest  flow  rate 
data  were  chosen  as  being  the  most  representative  of 
fully-developed  flow. 

Arc-column  electrical  characteristics  are  commonly 
presented  in  a  reduced  form  wherein  the  dependence  on  duct 

radius  (R)  is  removed  by  plotting  I/R  or  E  I  as  a  function 

Til  12  331 

of  E  RL   '   '    .   This  reduced  form  is  only  applicable, 

however,  under  those  conditions  for  which  both  non- 
equilibrium  effects  and  radiation  are  negligible.   Accord- 
ing to  the  variable  tube  radius  experiments  reported  by 
Emmons     and  Emmons  and  Land     and  according  to  com- 
parisons made  in  this  study  between  experiment  and 
equilibrium  and  nonequilibrium  numerical  solutions 
(Chapter  IV) ,  conditions  exist  for  which  radiation  and 
departures  from  equilibrium  have  a  significant  effect  on 
overall  arc  characteristics.   For  this  reason,  the  experi- 
mental results  which  were  selected  for  comparison  in 
Figure  III.  6  pertain  only  to  1.0  cm  diameter  ducts. 
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III. 3   Wall  Heat  Flux 

Typical  axial  distributions  of  the  total  wall  heat 
flux  (q)  for  several  arc  currents  and  laminar  gas  flow 
rates  in  argon  are  shown  in  Figures  III. 7  and  III. 8, 
respectively.   As  indicated  in  Figure  III. 8,  q  approaches 
a  constant  value  which  is  independent  of  gas  flow  rate. 
This  behavior  provides  additional  experimental  verifica- 
tion in  support  of  the  assumption  of  negligible  viscous 
dissipation. 

The  asymptotic  values  of  the  total  wall  heat  flux 

(q  )  for  Ar,  N2 ,  and  He  and  their  dependence  on  arc 

current  are  plotted  in  Figure  III. 9.   In  comparison  with 

the  values  for  argon,  the  larger  values  of  q   (and 

E    )  measured  in  nitrogen  and  helium  reflect  the  dis- 
x  /  a 

sociation  (N„)  and  relatively  high  ionization  (He) 

energies  which  characterize  these  two  gases.   Further 

comparisons  between  these  three  gases  are  deferred  to 

Chapter  IV. 

Because  of  the  absence  of  data  in  the  literature, 

comparison  of  the  total  wall  heat  flux  data  reported 

in  this  study  with  that  of  other  investigators  is  not 

possible;  however,  confidence  in  the  validity  of  the  q= 

a 

data  in  Figure  III. 9  may  be  established  from  the  favorable 

comparisons  of  Figure  III. 6  (Ev   versus  I).   An  energy 

x ,  a 
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balance  performed  on  a  segment  in  the  asymptotic  region 
dictates  that  the  rate  of  energy  removal  through  the 
tube  wall  must  be  balanced  by  the  local  input  rate  of 
electrical  energy.   That  is, 


q   •  2ttR  Ax  =  E     Ax  •  I .  (III.l) 

^a  x,a 


According  to  this  expression  the  asymptotic  value  of  the 
total  wall  heat  flux  is  expressed  by, 


<*a  =  T#  •  (III-2> 


Therefore,  confirmation  of  the  q  data  shown  in  Figure 
III. 9  is  implicit  with  the  favorable  comparisons  between 
the  asymptotic  electric  field  intensity  data  for  argon 
shown  in  Figure  III. 6.   Moreover,  the  local  energy  balance 
performed  on  a  segment  in  the  fully-developed  region  pro- 
vides a  convenient  check  on  the  accuracy  of  the  calormetric 
and  wall  potential  measurements.   As  is  evident  from  the 

q  and  E   I/2ttR  distributions  for  argon  shown  in  Figure 
a      x ,  a 

III. 9,  closure  of  this  local  energy  balance  is  obtained  to 
within  the  combined  experimental  uncertainties  of  q_(±  6%), 
Ev  a(±  3%)  ,  and  I(±  1%)  . 
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The  percentage  contribution  to  the  total  wall  heat 
flux  due  to  radiation  from  a  fully-developed  argon  plasma 

and  its  variation  with  arc  current  is  shown  in  Figure 

r  47  "i 

III. 10.   Also  shown  are  the  experimental  data  of  Barzelay 

[12  131 
and  Emmons    '    .   The  Harvard  data  are  actually  somewhat 

[13] 

redundant  in  that  the  distribution  reported  by  Emmons 

was  obtained  by  applying  a  correction  factor  of  1.5  to  the 
data  reported  by  Barzelay  (the  total  radiation  meter  used 
at  Harvard  is  discussed  in  the  Appendix) .   In  computing 
the  electrical  conductivity  and  radiation  source  strength 
properties  for  argon,  Emmons  determined  that  the  correction 
factor  was  necessary  in  order  to  obtain  the  best  fit  to 
overall  arc  data.   According  to  Figure  III. 10,  however, 
better  agreement  exists  between  the  radiation  data  of  this 
study  and  the  original  Harvard  data  reported  by  Barzelay. 
As  discussed  in  Chapter  II,  the  axial  distribution 
of  the  mixed-mean  enthalpy  (h  )  is  obtained  by  closing  off 
the  energy  balance  at  successive  axial  locations.   The 
corresponding  mean  temperatures  (T  )  are  then  obtained 
from  the  caloric  equation  of  state  [38,39,41].   Figures 
III. 11  and  III. 12  depict  the  respective  arc  current  and 
gas  flow  rate  dependency  of  the  axial  distribution  of  the 
mean  temperature  in  argon.   The  convergence  of  the  data 
in  Figure  III. 12  provides  further  evidence  of  the  attain- 
ment of  a  fully-developed  condition  which  is  independent 
of  operating  gas  flow  rate. 
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The  variations  with  arc  current  of  the  asymptotic 
mean  temperatures  in  Ar,  N2,  and  He  are  shown  in  Figure 
III. 13.   The  extraordinarily  low  value  of  the  mean  tempera- 
ture in  helium  (nearly  an  order  of  magnitude  lower  than 
the   corresponding  value  for  argon)  will  be  discussed  in 
Chapter  IV. 

The  Nusselt  number  (Nu  =  hD/k)  is  introduced  in  an 
attempt  to  correlate,  in  terms  of  appropriate  dimension- 
less  parameters,  the  wall  heat  transfer  in  the  asymptotic 
region  of  the  laminar,  constricted  tube  arc.   By  defining 
the  heat  transfer  conductance  (h)  in  terms  of  the  total 
wall  heat  flux  and  a  temperature  potential,  the  Nusselt 
number  takes  the  form 


2R  q 

Num  =  k  (T   -T  )  (III'3) 

m  m,a  w 


where  k   is  the  thermal  conductivity  evaluated  at  T   „ 
m  J  m  ,a 

and  T   is  the  duct  wall  temperature  (assumed  constant  at 
w 

500°K)*   Introduction  of  the  local  energy  balance  in  the 
asymptotic  region  (Equation  III. 2)  into  Equation  III. 3 
yields 


* 

While  this  value  of  T  was  chosen  so  as  to  be  consistent 

w 

with  the  assumed  value  used  in  the  equilibrium  solution 
by  Bower  and  Incropera  [20] ,  it  is  noted  that  a  consider- 
able variation  (say,  ±50%)  in  T  will  have  little  effect 

1  w 

on  the  difference  (T   -T  )  since  T    is  of  the  order  of 
1q4ok  ni,a  w         m,a 
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m  m,a  w 


The  electric  field  intensity  is  in  turn  related  to  the  arc 
current  by  an  expression  of  the  form 


R 

I  =  2ttE    J  ardr. 
x'a  o 


If  the  electrical  conductivity  (a)  is  assumed  to  be  con- 
stant and  is  evaluated  at  T    ,  then 

m,a 


I  *  ttR2E   a   .  (III. 5) 

x,a  m 


Substituting  this  relationship  into  Equation  III.  4  gives 
an  expression  of  the  form 


Nu   =   5 = .  (III.  6) 

(TTR)^a  k  (T   -T  ) 
mm  m, a  w 

It  is  now  convenient  to  introduce  an  Ohmic  heating  param- 
eter defined  as 


R2a  k  (T   -T  ) 
_,   _    mmm,aw  ,____. 

Oh,  =  ~— '- /  (III. 7) 

m         rz 
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which  may  then  be  used  to  correlate  the  Nusselt  number 
through  the  relationship 


Nu  =  ~  (Oh  )"1,0  .  (III.  8) 

m   „  2     m 

IT 


This  correlation  is  plotted  in  Figure  III. 14  for  argon 
along  wth  the  correlation  obtained  by  Bower  and  Incropera 
from  numerical  considerations.   While  the  numerically 
predicted  correlation  compares  reasonably  well  with  the 
experimental  data,  it  consistently  overpredicts  Nu  . 
Better  agreement  with  experiment  is  gained  by  multiplying 
Equation  III. 8  by  a  constant.   In  so  doing  the  experimental 
data  (for  argon  with  I  -  100  amp)  is  correlated  to  within 
4%  of  the  relationship 


Nu  =  1'\45    (oh  )"1,0  .  (III. 9) 

m     2      m 


Since  Equation  III. 8  is  independent  of  the  type  of  gas, 
the  above  correlation  should  apply  equally  well  for  nitro- 
gen and  helium.   While  the  experimental  values  of  the 
Nusselt  number  for  these  two  gases  (as  well  as  for  the 
lower  current  values  in  argon)  fall  within  the  range  of 

values  included  in  Figure  III. 14,  the  corresponding  values 

-3      -5 
of  Oh   fall  well  below  (-10    to  10   )  those  appearing  in 
m 

the  figure.   The  cause  of  these  rather  marked  deviations 
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from  Equation  III. 8  is  most  likely  due  to  the  inaccuracies 

which  result  from  evaluating  k  and  a   at  the  mean  temper- 

m      m 

ature  and/or  the  inaccuracies  which  result  from  use  of 
these  equilibrium  properties  under  conditions  which  are 
most  likely  characterized  by  significant  departures  from 
local  thermodynamic  equilibrium.   The  latter  point  will  be 
the  subject  of  further  discussion  in  the  following  chapter. 

III. 4   Anode  Heat  Flux 


Figure  III. 15  illustrates  the  dependency  of  the 
total  anode  heat  flux  on  arc  current  for  Ar ,  N2,  and  He. 
At  gas  flow  rate  and  arc  current  combinations  for  which 
the  anode  is  preceded  by  a  fully-developed  region,  the 
anode  heat  transfer  was  observed  to  vary  only  slightly 
(-  10%)  with  gas  flow  rate  and  to  vary  linearly  with  arc 
current  (for  I  >  10  0  amp)  according  to  the  expression 


Q„  =  8800  +  155  I  (III. 10) 

A 


for  argon ,  and 


QA  =  -  4300  +  367  I  (III. 11) 


for  nitrogen. 
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Figure  III. 15  Variation  of  Anode  Heat  Transfer  with 

Arc  Current 


70 


CHAPTER  IV   INTERPRETATION  OF  EXPERIMENTAL  RESULTS 
AND  COMPARISON  WITH  THEORETICAL  STUDIES 

IV. 1   Discussion 

The  performance  of  any  comparison  between  experi- 
mental and  theoretical  results,  of  course,  necessitates 
compatibility  between  the  experimental  configuration  and 
the  theoretical  model.   In  this  study,  this  compatibility 
is  achieved  only  for  the  asymptotic  region  of  the  arc. 
Although  considerable  data  have  been  acquired  for  the  in- 
let region  of  the  constrictor,  it  was  determined  that  mea- 
surements made  in  this  region  were  extremely  sensitive  to 
the  entrance  configuration.   In  particular,  variation  of 
the  duct  inlet  geometry  (effected  through  adjustment  of 
the  cathode-nozzle  gap)  significantly  influenced  the  mea- 
sured axial  distributions  of  the  electric  field  intensity 
and  total  wall  heat  flux  in  the  inlet  region.   Similarly, 
all  theoretical  considerations  for  the  arc-heating  region 
require  the  assumption  of  entrance  flow  profiles  and,  as 
shown  by  Bower  and  Incropera     ,  predictions  for  the  in- 
let region  depend  upon  the  assumed  form  of  the  initial 
velocity  and  enthalpy  profiles.   Accordingly,  this  sensi- 
tivity of  experimental  and  theoretical  results  to  entrance 
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conditions  precludes  any  meaningful  comparison  for  the  in- 
let region.   In  contrast,  the  experiments  of  this  study 
and  the  numerical  solution     indicate  that  entrance  ef- 
fects have  no  bearing  on  either  measured  or  predicted  quan- 
tities for  the  asymptotic  region  and,  therefore,  a  meaning- 
ful comparison  of  results  applicable  to  this  region  is 
possible. 

This  chapter  will  deal  only  with  comparisons  between 
experiment  and  theory  which  apply  to  the  asymptotic  (and 
anode)  regions  of  the  constricted  tube  arc.   In  particular, 
comparisons  for  argon  will  be  made  between  the  experimental 
results  of  this  study  and  the  analytical  solution  by  Stine 

ri8i 

and  Watson     ,  the  numerical  equilibrium  solution  by 

Bower  and  Incr opera     ,  and  the  numerical  nonequilibrium 

T251 
solution  by  Clark     .   Owing  to  the  scarcity  of  experi- 
mental data  acquired  in  this  study  for  helium,  the  experi- 

F131 
mental  data  of  Emmons     will  be  used  for  comparison  with 

the  theoretical  results  obtained  by  Bower     ,  Pytte  and 
Williams     ,  and  Pytte  and  Winsor     .   For  nitrogen,  com- 
parisons will  be  made  between  the  experimental  results  of 
this  study  and  the  Stine-Watson  analytical  solution.   In 
making  these  comparisons,  emphasis  will  be  placed  upon 
ascertaining  the  influence  of  nonequilibrium  effects  on  the 
arc  parameters.   In  comparing  the  anode  heat  transfer  model 

r  6  7i 

of  Pfender,  et  al l  '  J  with  the  experimental  results  of 
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this  study,  attention  is  given  to  determining  the  ability 
of  this  model  to  accurately  predict  the  anode  heat  flux, 
particularly  in  light  of  available  empirical  correlations. 

IV. 2   Comparisons  For  Argon 
Although  the  equilibrium  model  of  Bower  and  Incro- 
pera  and  the  nonequilibrium  model  of  Clark  have  been  ad- 
equately discussed  previously,  more  must  be  said  of  the 

Stine-Watson  analytical  solution  before  using  it  in  sub- 

[5] 
sequent  comparisons.   As  reported  by  Jedlicka  and  Stine    , 

the  expressions  for  the  electric  field  intensity  and  total 

wall  heat  flux  which  result  from  this  closed  form  solution 

are  given  by 


C,     -C9L/m  ~1/2 
Ex  =  g±  [1-e  4         ]  (IV. 1) 


and 


_       -C.L/m  1/2 
q  =  C,  ~  [1-e   Z    ]     ,  (IV. 2) 

*  d^ 


where  the  constants  C2,  C3,  and  C.  are  specific  functions 
of  the  constant  pressure  gas  transport  properties,  d  is 
the  duct  diameter,  and  L  is  the  duct  inlet  length.   When 
applied  to  the  asymptotic  region  (L  ■*  °°)  of  the  constrict- 
ed tube  arc,  Equations  IV. 1  and  IV. 2  reduce  to  the  expres- 
sions : 
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Ex,a  =  IT  <IV-3» 


9a  "  C4  72   •  (IV-4) 

d 


Basing  the  above  equations  upon  the  tube  radius  and  using 

[49] 
the  following  results  for  the  constants  C3  and  C. 


CT  =  ^T77  (IV. 5) 

J    hg   ' 


C4  =  °4*§  .  (IV. 6) 


Ag" 


Equations    IV. 3    and   IV. 4    take    the    form 


Va  =  ^T72i  <IV-7> 


and 


qa  =  — T72  ~2  (IV. 8) 

a   Ag1/2  R2 


Ag  (mho/watt)  represents  the  assumed  linearized  relation- 
ship between  the  electrical  conductivity  and  the  heat  flux 

potential  <f>  (<J>  =  JkdT)  ,  and  the  units  for  the  remaining 

2 

variables  are  E    (volts/m) ,  q  (watts/m  ),  I (amp),  and 
x  /  a  a 
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r  i  o  i 
R(m) .   While  the  results  reported  by  Stine  and  Watson 

are  primarily  for  nitrogen  and  air,  this  simplified  model 

has  been  modified  by  Bower     to  include  the  appropriate 

linearized  property  relationships  for  argon.   By  applying 

the  value  of  Ag  reported  by  Bower,  the  simplified  Stine- 

Watson  expressions  (for  argon  with  R  =  0.5  cm)  are  then 

given  by 


E     =  9.55  volts/in  (IV. 9) 

x,a 


and 


q   (Btu/ft2-min)  =63.4  I (amp)  .         (IV. 10) 
a 


The  asymptotic  values  of  the  electric  field  inten- 
sity (E    ) ,  total  wall  heat  flux  (q  ) ,  fraction  of  radia- 
x ,  a  a 

tive  to  total  wall  heat  flux  (q   ,   /q  ) ,  mixed-mean  en- 

raci ,  a   a 

thalpy  (h    ) ,  and  mean  temperature  (T    ) ,  as  computed 
m ,  a  m ,  a 

from  the  equilibrium     and  nonequilibrium     models, 
are  listed  in  Table  IV. 1  along  with  the  corresponding  ex- 
perimental results  for  an  argon  arc.   The  values  from 
Table  IV. 1,  along  with  results  obtained  from  Equations 
IV. 9  and  IV. 10  for  the  Stine-Watson  model,  were  used  to 
generate  the  plots  shown  in  Figures  IV. 1  through  IV. 5. 

The  large  disparities  (in  both  the  magnitude  and 
the  functional  dependence  on  arc  current)  between  the 
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E     (Figure  IV.  1)  and  cr   (Figure  IV. 2)  distributions  com- 

X  /  3.  a. 

puted  from  the  simplified  Stine-Watson  model  and  those  ob- 
tained experimentally  are  thought  to  be  the  result  of  two 
oversimplifications  inherent  in  the  model.   In  particular , 
the  use  of  linearized  property  relationships  and  the  as- 
sumption of  negligible  radiation  losses,  particularly  at 
the  higher  arc  currents ,  are  thought  to  be  the  primary 
cause  of  these  disparities. 

The  comparison  between  the  numerical  solution  by 
Bower  and  Incropera  and  the  experimental  results  indicate 
that  the  equilibrium  theory  agrees  favorably  with  experi- 
ment in  the  100  amp  range,  while  it  overpredicts  and  under- 

predicts  the  experimental  data  for  both  E    and  q   at 

x  ,  a      a 

lower  and  higher  arc  currents,  respectively.   A  plausible 
explanation  for  the  overprediction  at  the  lower  arc  cur- 
rents can  be  obtained  by  considering  the  following  expres- 
sion for  Ohm's  law 

E   =  -gS .  (IV. 11) 

x   2tt  /k  ardr 
o 

If  the  low  current  arc  is  characterized  by  significant  de- 
partures from  thermochemical  equilibrium,  the  electrical 
conductivity  used  in  the  equilibrium  solution  will  be  less 
than  the  actual  value  of  a   in  those  regions  where  nonequi- 
librium  conditions  prevail.   Consequently,  the  theoretical 
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Figure  IV. 1  Comparison  Between  Experiment  and  Theory 
for  the  Asymptotic  Electric  Field  Inten- 
sity in  Argon 


79 


<D 

H 

to 

-P 
O 

Ei 

U 

OJ 

•H 

ro 

■P 

S 

00 

< 

OJ 

(1) 

M 
O 

<fr 

IH 

cvj 

P. 

0 

to 

i 

(1) 

O 

Eh 

o 

rvi 

X 

'd 

*""■* 

C 

C 

f0 

E 
i 

<x> 

CM 

i 

\ 

•H 

3 

00 

a.  o 

OJ 

o 

Between 
Flux   in 

00 

c 

0  -P 

•h  a) 

H  SB 
£3  -h 

«■ 

(dujD)    luejjno    ojv 


(N 


u 

•H 


80 


evaluation  of  the  integral  of  the  electrical  conductivity 
over  the  tube  cross-section  will  be  reduced,  which,  ac- 
cording to  Equation  IV. 11,  will  result  in  an  overprediction 

of  E    ,  increased  Joule  heating,  and  thus  an  accompanying 
x ,  a 

overprediciton  of  the  total  wall  heat  flux.   While  this 

nonequilibrium  effect  suitably  explains  the  lower  current 

comparisons,  it  does  not  explain  why  the  equilibrium  theory 

apparently  underpredicts  E    and  q   at  the  higher  arc 

x ,  a      a 

currents . 

Recent  theoretical  studies  by  Kruger     ,  Incropera 

r  23 1  T2R1 

and  Viegas     ,  and  Clark     and  recent  experimental 

[211 
spectroscopic  data  reported  by  Bott     and  by  Incropera 

r  221 

and  Giannaris      indicate  that  arc  nonequilibrium  effects 

are  most  severe  at  low  arc  currents  and  diminish  with  in- 
creasing arc  current.   On  this  basis,  the  experimental  re- 
sults and  equilibrium  predictions  which  are  compared  in 
Figures  IV. 1  and  IV. 2  should  converge  with  increasing  arc 
current  when,  in  fact,  they  diverge  for  I  >  100  amp.   Ac- 
cordingly, the  equilibrium  underpredictions  for  E    and 

x ,  a 

q  at  higher  arc  currents  would  not  appear  to  be  a  result 
of  wall-induced  nonequilibrium. 

A  plausible  explanation  for  the  high  current  dis- 
parity is  the  inability  of  the  equilibrium  model  to  prop- 
erly account  for  conduction  effects  at  the  elevated  temp- 
eratures associated  with  high  current  operation.   In 
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particular,  it  is  possible  that  the  equilibrium  model  is 
underpredicting  the  energy  transfer  to  the  tube  wall  due 
to  conduction.   This  contention  is  supported  by  the  com- 
parisons for  the  fraction  of  radiative  to  total  wall  heat 
flux  shown  in  Figure  IV. 3,  as  well  as  by  the  mean  enthalpy 
and  mean  temperature  comparisons  of  Table  IV. 1.   Specif- 
ically, it  is  noted  from  Figure  IV. 3  that,  while  excellent 
agreement  exists  between  experiment  and  the  nonequilibrium 

predictions  for  q   ,   /q  ,  the  equilibrium  prediction,  for 

rao. ,  a   a 

I  >  100  amp,  is  nearly  double  the  nonequilibrium  and  ex- 
perimental distributions.   In  comparing  the  equilibrium 
and  nonequilibrium  predictions  for  q   ,   /q  ,  it  may  be 

1C  did  f  cl         3. 

argued  that  the  observed  differences  are  possibly  a  re- 
sult of  a  dissimilarity  in  the  radiation  source  density 

terms  used  in  the  two  solutions.   This  is  not  the  case, 

[251 
however,  for  it  has  been  shown  by  Clark     that,  in  the 

limit  of  equilibrium,  the  radiation  model  employed  in  the 
nonequilibrium  solution  agrees  favorably  with  the  radia- 
tion source  density  used  by  Bower  and  Incropera.   More- 
over, both  of  these  plasma  radiation  models  are  consistent 
with  the  assumed  wavelength  distribution  of  the  radiation 
incident  upon  the  surfaces  of  the  radiation  gage  (see  the 

Appendix) .   The  overprediciton  of  q   ,   /q   is  therefore 

rao. ,  a   a 

most  likely  due  to  the  higher  temperature  levels  which 
characterize  the  equilibrium  solution.   This  fact  is 
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further  evidenced  by  the  T    comparisons  of  Table  IV. 1  as 

m,a    c 

well  as  by  the  comparisons  between  equilibrium,  nonequilib- 

rium,  and  experimentally  determined  temperature  pro- 

*.,   [25] 
files 

On  the  basis  of  these  arguments,  it  is  concluded 

that,  for  the  higher  (I  >  100  amp)  arc  currents  in  argon, 

the  reduction  of  the  equilibrium  temperature  profiles 

needed  to  bring  about  a  corresponding  reduction  in  the 

electrical  conductivity  and  thus  an  increase  in  E    and 

x  x,a 

q   (along  with  a  decrease  in  q   ,   /q   and  T    )  could  be 
a  raQ,a   a      m,a 

achieved  if  the  equilibrium  model  were  to  provide  for 
greater  energy  loss  by  conduction.   In  an  effort  to  con- 
firm this  premise,  some  attention  was  given  to  the  thermal 

r  3  7 1 

conductivity  model  of  DeVoto     which  was  employed  in  the 

equilibrium  solution. 

In  addition  to  the  conventional  translational  mode 
of  energy  conduction,  there  will  be  a  transport  of  chem- 
ical energy  (in  the  form  of  the  energy  of  recombination) 
through  the  ordinary  (or  concentration)  diffusion  of  elec- 
tron-ion pairs  towards  the  outer  regions  of  the  arc. 
While  the  global  conservation  equations  employed  in  the 
equilibrium  model  do  not  explicitly  account  for  ordinary 
diffusion,  an  attempt  is  made  to  implicitly  account  for 
the  resulting  transport  of  ionization  (or  recombination) 
energy  through  use  of  a  reactive  thermal  conductivity 
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r  37 1 

(X  ) .   From  the  work  of  DeVoto    *,    X_  is  defined  by  the 
heat  flux  vector  expression 


q  =  -k'VT  -  X_VT  ,  (IV. 12) 

where  k1  is  the  translational  thermal  conductivity.   This 
expression  evolves  from  the  more  general  form  of  the  heat 
flux  vector  equation  by  replacing  the  diffusion  contribu- 
tion to  q  by  the  term  -A  VT,  wherein  it  is  assumed  that 
the  gradient  of  concentration  can  be  replaced  by  a  temper- 
ature gradient  through  the  equilibrium  law  of  mass  action*. 
The  point  to  be  made  here  is  that  the  critical  assumption 
inherent  in  DeVoto 's  formulation  of  the  reactive  thermal 
conductivity  is  that  of  local  chemical  equilibrium.   More 
specifically,  it  is  assumed  that,  in  the  case  of  a  chem- 
ically reacting  mixture,  all  forward  and  reverse  reaction 
rates  are  sufficiently  fast  to  preclude  any  effect  on  the 
local  composition  due  to  ordinary  diffusion.   Otherwise, 
use  of  the  equilibrium  law  of  mass  action  to  relate  con- 
centration and  temperature  gradients  could  not  be  justified. 
Conversely,  in  a  constricted  tube  plasma  which  is  charac- 
terized  by  strong  concentration  gradients,  it  is  certainly 


This  derivation  is  treated  in  detail  by  Butler  and 
Brokaw[51]  and  by  Fayt52]. 


85 


possible  that  the  local  atomic  composition  will  be  signif- 

r  25 1 

icantly  influenced  by  ordinary  diffusion     ,  in  which 
case  the  assumption  of  local  chemical  equilibrium  is  in- 
valid.  It  is  therefore  suggested  that  the  high  current 
disparity  between  experimental  and  equilibrium  theoretical 
results  is  due,  in  part,  to  a  chemical  nonequilibrium  ef- 
fect on  the  reactive  thermal  conductivity. 

In  connection  with  the  above  discussion,  it  is  im- 
portant to  note  that  this  observation  is  seemingly  incon- 
sistent with  previous  considerations  which  indicated  that 
the  departure  from  thermochemical  equilibrium  is  less  se- 
vere at  higher  arc  currents.   However,  in  comparing  DeVoto ' s 

calculations  for  the  translational  and  reactive  thermal 

T371 
conductivities     ,  it  is  noted  that  the  percentage  contri- 
bution to  the  effective  thermal  conductivity  (kEFF  =  k'+A  ) 
due  to  A   increases  from  zero  at  T  =  6,000°K  to  approxi- 
mately  50%  at  T  =  14,000°K.   Accordingly,  any  errors  which 
may  be  introduced  as  a  result  of  using  DeVoto ' s  formulation 
for  A   would  become  increasingly  more  significant  at  higher 
arc  temperatures.   In  this  respect,  it  is  suggested  that  the 
failure  of  the  equilibrium  solution  to  provide  for  suffi- 
cient energy  loss  by  conduction  at  the  higher  arc  currents 
(I  >  100  amp)  is  due  to  both  its  inability  to  allow  for 
sufficient  transport  of  ionization  energy  to  the  tube  walls 
and  its  inability  to  explicitly  account  for  the  effects  of 
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diffusion  on  the  local  gas  composition.   In  particular,  it 
is  suggested  that  the  thermal  conductivity  values  (specif- 
ically, X    )    employed  in  the  equilibrium  solution  should 
be  higher  in  order  to  bring  the  equilibrium  theory  predic- 
tions into  better  agreement  with  both  the  experimental 

observations  and  the  nonequilibrium  predictions.   This  con- 

[15] 
tention  is  supported  by  the  statements  of  Maecker     and 

T141 
by  the  comparisons  reported  by  Morris,  et  al.     .   These 

indicate  a  significant  increase  (by  as  much  as  a  factor 
of  two)  in  experimentally  determined  values  of  thermal 
conductivity  over  the  corresponding  predictions  by  DeVoto 
for  temperatures  greater  than  10,000°K. 

While  attempts  have  been  made  to  interpret  Figures 
IV. 1  through  IV. 3  on  the  basis  of  nonequilibrium  effects 
and  thermal  conductivity  considerations,  final  judgement 
on  the  adequacy  of  the  equilibrium  model  cannot  be  made 
until  the  combined  experimental  and  theoretical  uncertain- 
ties are  considered.   In  an  attempt  to  determine  the 

range  of  uncertainties  associated  with  the  equilibrium 

[53] 
solution,  a  parametric  study  was  conducted     wherein 

the  properties  (a,  kEFF/  and  P  )  were  selectively  varied 

over  their  respective  limits  of  estimated  accuracy 

(a,  ±  15%;  k„__,  ±  10%;  and  PD,  ±  15%)*.   It  was  determined 


*The  theoretical  uncertainties  for  a  and  k___  were  set  by 

[60] 
DeVoto  (private  communication  with  Kingsbury    )  while 

the  uncertainty  in  P  was  arbitrarily  selected 

R 
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from  this  study  that  the  upper  uncertainty  limit  for  E 

x ,  a 

and  q=  resulted  from  decreasing  a  by  15%,  increasing  k„__ 

a  h,r  r 

by  10%,  and  increasing  P   by  15%,  while  the  lower  limit 

resulted  from  increasing  a  by  15%,  decreasing  k„„„   by  10%, 

tit  t 

and  decreasing  P  by  15%.   The  corresponding  uncertainty 

bands  for  the  equilibrium  predictions  are  shown  in  Figures 

IV. 4  and  IV. 5  along  with  the  experimental  data  and  the 

related  experimental  uncertainities ,  E     (±3%)  and 

x ,  a 

q  (±6%) .   On  the  basis  of  the  overlap  between  the  theoret- 
ical and  experimental  uncertainty  bands ,  and  in  view  of 
the  previous  discussion  concerning  the  use  of  DeVoto ' s 
formulation  for  the  thermal  conductivity,  it  is  concluded 
that,  to  within  the  existing  knowledge  of  plasma  transport 
and  radiation  properties,  single-fluid  models  are  suitable 
for  predicting  asymptotic  values  for  the  electric  field 
intensity  and  total  wall  heat  flux,  at  least  for  the  range 
of  conditions  in  argon  considered  in  this  study.   Converse- 
ly, the  equilibrium  model  is  not  suitable  for  predicting 

q   ,   /q  over  the  entire  range  of  conditions  investigated, 
r  ao ,  a   a 

nor  is  it  adequate  for  predicting  h    and  T   „  at  the 
^         c        3   m,a      m,a 

lower  (<  50  amp)  and  higher  (>  150  amp)  arc  current  con- 
ditions.  In  addition,  it  is  felt  that  the  equilibrium 

model  would  not  be  suitable  for  predicting  E^    and  q 

x ,  a      a 

if  the  above  comparisons  were  extended  to  include  condi- 
tions of  lower  and  higher  arc  currents.   In  this  respect, 
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Figure  IV. 4   Comparison  Between  the  Combined  Experimen- 
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it  is  concluded  that  a  more  rigorous  theoretical  treat- 

r  25 1 

ment,  such  as  the  nonequilibrium  solution. by  Clark     , 

would  be  required  to  accurately  model  the  asymptotic  re- 
gion of  the  constricted  tube  argon  arc . 

Before  concluding  this  section  on  comparisons  for 
the  asymptotic  region  in  argon,  some  explanation  is  in 
order  concerning  the  lack  of  agreement  between  the  exper- 
imental value  and  the  nonequilibrium  prediction  for  h 

at  I  =  35  amp  (Table  IV. 1) .   Calculation  of  h    in  the 

m  ,  a 

r  25 1 

nonequilibrium  solution     is  effected  through  the  expres- 

•   [54] 
sion 


JR  puhrdr 

K   a  E  "5 '  (IV. 13) 

m,a    rR    , 
'     J   purdr 

where  p(r)  is  the  local  density  of  the  plasma  and  u(r)  is 
the  global  axial  velocity  component.  The  radial  enthalpy 
distribution,  h(r),  required  for  the  above  calculation  is 
given  by 

I. 

h(r)  =  §  g-  [aT  (r)  +  Th(r)]  +  -JpH  a(r)  ,    (IV. 14) 
A  A 

where  the  first  two  terms  represent  the  contributions  due 
to  the  electrons  and  heavy  particles,  respectively,  and 
the  remaining  term  represents  the  chemical  enthalpy  of 
recombination  (the  specific  variables  are  defined  in  the 
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nomenclature).   In  formulating  Equation  IV. 14,  the  addi- 
tional enthalpy  contributions  due  to  internal  electronic 
excitation  and  radiation  source  density  effects  were  as- 
sumed negligible .   The  governing  equations  for  the  asymp- 
totic region  are  obtained  from  the  general  set  of  conser- 
vation equations  by  setting  all  axial  derivatives  (with 

the  exception  of  the  pressure  gradient,  dp/dx)  equal  to 

[251 
zero     .   This  manipulation,  along  with  the  assumptions 

of  negligible  viscous  dissipation  and  pressure  work,  re- 
sults in  a  decoupling  of  the  hydrodynamic  and  thermal  flow 
fields.   Consequently,  the  velocity  profiles  required  in 

the  above  definition  for  h    can  only  be  obtained  by 

m,a 

solving  the  following  asymptotic  axial  momentum  equation 


I  |_  [yr  |Uj  _  d£  =  (IV. 15) 

r  8r      8r    dx 


It  is  the  uncertainty  associated  with  the  determi- 
nation of  dp/dx  in  Equation  IV. 15  that  renders  the  non- 
equilibrium  solution  unable  to  accurately  predict  u(r) 

and  thus  h    at  arc  currents  less  than  50  amp.   For  the 
m,a 

same  reason,  nonequilibrium  predictions  for  h    at 
currents  greater  than  100  amp  are  unrealistic  -  the  agree- 
ment in  the  50  to  100  amp  range  is  thought  to  be  fortu- 
itous . 
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IV. 3   Comparisons  for  Helium  and  Nitrogen 
Asymptotic  values  for  helium  of  the  electric  field 

intensity,  total  wall  heat  flux,  mean  enthalpy,  and  mean 

[551 
temperature  computed  from  the  equilibrium  model  by  Bower 

are  listed  in  Table  IV. 2.   Also  listed  in  this  table  are 

the  experimental  data  for  helium  reported  by  McKee ,  et 

r  3  2 1  T131 

al.     and  Emmons     ,  as  well  as  the  50  ampere  data  of 

this  study.   The  data  from  Table  IV. 2  were  used  for  the 

plots  of  Figures  IV. 6  and  IV. 7,  wherein  it  is  noted  that 

the  experimental  q  values  (with  the  exception  of  the 

a 

single  q  value  measured  calorimetrically  in  this  study) 
a 

were  obtained  from  the  asymptotic  energy  balance  equation 
(Equation  III. 2)  and  the  experimental  E    data. 

X  f  Si 

In  comparing  the  equilibrium  theory  predictions 
with  the  experimental  results,  it  is  readily  apparent 
that,  unlike  the  previous  comparisons  for  argon,  the 
single-fluid  model  is  totally  deficient  in  its  attempt 
to  accurately  describe  the  overall  energetics  of  the  heli- 
um arc.   The  reason  for  this  deficiency  is  believed  to  be 
a  result  of  the  severe  nonequilibrium  conditions  which 
are  known  to  exist  in  low  current  (I  <  400  amp) ,  atmo- 
spheric helium  arcs. 

The  pronounced  departure  of  helium  arc  plasmas 
from  conditions  of  local  thermodynamic  eqiulibrium  has 
been  observed  and  discussed  by  others,  among  the  first  of 


93 


-OS 

go 
En     — 


e 

(d  .H 
•»\ 

X!      -P 

pq 


-p 

^ 

*      iw 

C 

fd\ 

•H 

O1     3 

6 

-P 

PQ 

l 

03 

fd  -P 

-iH 

c 

X   0 

-H 

W      > 

CD 
O 

o 

CO 


a 


o 
o 

I    00 

I     •* 
ro 


o 

CM 


o  o 

O  IT) 

in  oo 


rr  <j\ 


^r  ^ 


CM  (N 


in 
in 


CM 

ro 


4-1  nH 

c  o 

CD  CO 
•H     • 

CD  -H 

a  3 


in 
ro 


o  o 
in  o 
ro    I   r» 


i 


in 


o       o 
co       o 


i 


ON 


o 

o 

I   m 

I     * 


o 
o 

I    CO 

I     * 

CM 


O  O  O 

o  o  o 

00  wo 


r*-  <y>  (N 


O    H    nH 


O  CT>  <£> 
vo  in  co 


in 

in 


ro 


•P 

c 

CD 

e 

•H 
>H 

<L) 


-P  H 
C  0 
0)  CO 

6 

•H  • 
J-l  -H 
CD  -H 

www 


o 
m 


o 

o 

i  vo 

i     » 

CM 


O 

O 

I    r- 

I  - 
iH 
ro 


O  O 

o  o 

CM  CM 


CM  <T» 
ro  iH 


•         • 

00  00 

^tf  rsi 


in 
in 


ro 


-P  H 

C  0 

CD  W 

£ 

•H  • 
U  -H 
CD  -H 

a  3 
x  u1 
w  w 


o 
o 


o  o 
o  o 
in  co 


ro  m 

^  CM 


co  r- 


ro  in 

^J<   CM 


in 
m 


ro 


P    rH 

c  o 

CD  CO 
6 
•H     • 

U    rH 

CD  -H 

a  3 

X    D1 
W  W 


o 

in 


o 
o 
in 


o 
o 


ro 
in 


ro 


o 


oo 
ro 


ro 


ro 


-P 
C 

CD 

6 

•H 
^1 
CD 
Oi 
X 

W 


o 

o 

CM 


■P 
C 

CD 

S 

•H 
U 

CD 

a 

X 

w 


o 
in 

CM 


S 

o 
•a 

CD 
P 

6 
0 
0 

CD 
U 
CD 
5 


ro 


C 
rd 


CM 

ro 


W 
CD 
U 
C 
CD 
U 
CD 
4-1 
CD 
« 

U 
0 

14-1 

W 
CD 

P 

rH 

<d 

> 

X 

a 

rH 
M-4 

P 

<d 

CD 


rd 
2 


rd 
p 
0 

9n 


94 


which  were  investigators  at  Harvard    '   '   '    .   These 
studies  were  primarily  concerned  with  determination  of 
the  plasma  transport  properties,  a  and  k,  through  use  of 

the  so-called  inversion  method.   As  discussed  by  Emmons 

[121 
and  Land     ,  this  method,  when  formulated  under  the  as- 
sumptions of  local  thermodynamic  eqiulibrium  and  negli- 
gible radiative  losses,  dictates  that  the  reduced  arc- 
column  characteristics  (E   I  versus  E   R)  should  be 

x,a  x,a 

related  by  a  universal  curve  independent  of  tube  radius. 

However,  in  performing  electric  field  intensity  measure- 

[13] 
ments  on  helium  arcs  of  varying  tube  radius ,  Emmons 

determined  that  the  parameters  E   I  and  E   R  were  not 

r  x,a       x,a 

uniquely  related.   Since  the  assumption  of  negligible 

radiation  is  apparently  valid  for  low  current  (I  <  400 

r  57 1 

amp),  atmospheric  helium  arcs     ,  Emmons  concluded  that 

helium  arcs  must  be  characterized  by  a  significant  de- 
parture from  conditions  of  local  thermodynamic  equilibrium. 
Further  evidence  of  severe  nonequilibrium  conditions 

in  helium  arcs  is  provided  by  Kruger     ,  Bott     ,  and 

T  58 1 
Aleksandrov,  et  al. L    .   Specifically,  these  investigators 

concluded  that  the  free  electrons  and  upper  excited 

atomic  states,  while  in  equilibrium  with  one  another,  are 

not  in  equilibrium  with  the  lower  excited  (n  <  3)  and 

ground  states.   In  addition,  Kruger  notes  that,  in  the 

outer  regions  of  the  arc,  the  electron  number  densities 
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may  be  several  orders  of  magnitude  higher  than  would  be 
predicted  if  the  local  composition  were  determined  from 
the  assumption  of  Saha  equilibrium.   The  existence  of 
significant  departures  from  equilibrium  can  also  be  in- 
ferred from  the  helium  data  acquired  in  this  study.   As 
indicated  in  Table  IV. 2,  energy  balance  measurements  per- 
formed on  the  helium  arc  dictate  an  asymptotic  mean  tem- 
perature of  1350°R  (or  750°K)  and  yet,  even  at  such  a 
comparatively  low  value  of  T    ,  there  is  sufficient  ion- 
ization to  produce  a  luminescent  arc  which  is  capable  of 
sustaining  a  current  of  50  amperes.   This  observation 
suggests  the  existence  of  a  severe  thermochemical  non- 
equilibrium  condition  in  the  helium  arc.   Similar  ob- 

r  58 1 

servations  were  made  by  Aleksandrov,  et  al .      whose 
spectroscopically  determined  electron  temperatures  were 
nearly  1.5  orders  of  magnitude  greater  than  their  mean 
temperatures  based  upon  energy  balance  measurements. 

On  the  basis  of  these  findings,  the  comparisons 

[55] 
between  experiment  and  the  equilibrium  theory  of  Bower 

(Figures  IV. 6  and  IV. 7)  are  virtually  meaningless.   That 
is,  since  the  plasma  tube  flow  situation  which  the  single- 
fluid  model  is  attempting  to  describe  bears  little  re- 
semblance to  the  severe  nonequilibrium  state  which  actu- 
ally exists,  a  systematic  discussion  of  the  differences 
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is  not  possible.   As  indicated  in  Figure  IV. 6,  the  pre- 
diction of  the  arc-column  characteristics  in  helium  by 
Pytte  and  Williams     offers  some  improvement  over  the 
numerical  solution  by  Bower;  however,  this  theoretical 

model  continues  to  underpredict  E    by  some  20%.   This 

x  ,  a 

underprediction  by  the  Pytte-Williams  model  is  most  likely 

f  391 

due  to  the  use  of  equilibrium  calculations  for  a  and  k l 

as  well  as  to  use  of  the  Elenbaas-Heller  equation.   Pytte 

[241 
and  Winsor     have  given  greater  consideration  to  non- 
equilibrium  in  atmospheric  helium  arcs  by  proposing  a  two- 
temperature  theory  which  is  based  upon  the  asymptotic 
moment  equations  of  continuity,  momentum,  and  energy  as 
derived  from  an  electron  transport  (Boltzmann)  equation. 
While  this  model  accounts  for  thermal  nonequilibrium  be- 
tween the  electrons  and  heavy  particles,  it  suffers  from 

[39] 
use  of  equilibrium  thermal  conductivity  values     and 

use  of  the  Saha  equation  (evaluated  at  the  local  electron 
temperature)  to  determine  the  local  composition.   In  ob- 
serving the  excellent  agreement  which  exists  between  the 
Pytte-Winsor  predictions  and  Emmons  *  experimental  data  for 

E     (Figure  IV. 6) ,  it  would  appear  that  use  of  the  Saha 
x ,  a 

equation  and  eqiulibrium  thermal  conductivities  is  justi- 
fied. This  agreement  is  misleading,  however,  for  Pytte  and 
Winsor  adjusted  their  electrical  conductivity,  a,  in  order 
to  obtain  the  best  fit  to  Emmons1  experimental  reduced 
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arc-column  data  for  p  =  1.0  atm  and  R  =  0.5  cm.   In  apply- 
ing the  resulting  value  of  a  to  calculations  for  tube 
radii  other  than  R  =  0.5  cm,  the  Pytte-Winsor  model  does 
predict  the  trends  indicated  by  the  experimentally  de- 
termined shifts  in  the  E    R  versus  E    I  curves;  however, 

x,a  x,a 

[24] 
it  fails  to  predict  the  magnitude  of  these  shifts     .   It 

is  believed  that  this  deficiency  results  from  the  failure 
of  the  Pytte-Winsor  model  to  properly  account  for  chemical 
nonequilibrium.   Accordingly,  the  need  remains  for  a  rig- 
orous, theoretical,  multifluid  model  which  treats  the 
thermochemical  nonequilibrium  flow  of  helium  in  the  inlet 
and  asymptotic  regions  of  the  constricted  tube  arc. 
Similarly,  the-  author  knows  of  no  multifluid 
analyses  in  the  literature  from  which  predicitons  of  the 
global  parameters,  E   and  q  can  be  made  for  the  inlet 
and  asymptotic  regions  of  nitrogen  arcs.   Consequently, 

the  following  comparisons  for  nitrogen  are  limited  to  the 

[  18  ] 
single-fluid  analytical  model  by  Stine  and  Watson 

As  was  the  case  for  argon,  the  generalized  Stine-Watson 

expressions  for  E    and  q   (Equations  IV. 7  and  IV. 8) 

x ,  a      a 

must  be  adapted  to  the  particular  gas  and  tube  radius 
being  considered.   For  this  purpose,  the  values  of  E 

X  f  d 

[3] 
and  R  for  nitrogen  reported  by  Shepard  and  Watson    were 

used  with  Equation  IV. 7  to  determine  Ag.   This  value  was 
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then  used  to  give  the  following  simplified  expressions  for 
the  Stine-Watson  model,  as  applied  to  the  asymptotic  re- 
gion of  a  nitrogen  arc  with  R  =  0.5  cm. 


Ev   =  26.4  volts/in  (IV. 16) 


qa  =  175.5  I  (IV. 17) 

These  relationships,  along  with  the  experimental 
data  from  Figures  III. 6  and  III. 9,  were  used  for  the  com- 
parisons of  Figures  IV. 8  and  IV. 9.   Relative  to  previous 
comparisons  for  argon,  application  of  the  Stine-Watson 
model  to  nitrogen  provides  better  agreement  with  exper- 
iment.  While  the  temperature  dependence  of  the  radiation 

source  densities  (P  )  for  argon  and  nitrogen  are  nearly 

r  59  "I 

the  same     ,  the  mean  temperatures  in  nitrogen  are  some 

50%  lower  than  the  corresponding  values  of  T    in  argon 

c  3  m,a 

(Figure  III. 13).   Therefore,  the  percentage  contribution 
to  the  total  wall  heat  flux  due  to  radiation  from  a  nitro- 
gen arc  plasma  would  be  considerably  less  than  the  corre- 
sponding contribution  for  an  argon  plasma  of  the  same  arc 
current.   As  such,  the  negligible  radiation  assumption 
implicit  in  the  Stine-Watson  model  is  less  severe  for 
nitrogen;  nevertheless,  as  shown  in  Figures  IV. 8  and  IV. 9, 
this  assumption  continues  to  render  this  model  inaccurate 
at  higher  arc  currents. 
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Figure  IV. 8  Comparison  Between  Experiment  and  Theory 
for  the  Asymptotic  Electric  Field  Inten- 
sity in  Nitrogen 
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In  making  further  comparisons  between  the  three 

gases,  it  is  noted  from  Figure  III.  6  that,  for  a  given  arc 

current,  successively  higher  voltages  are  required  to 

sustain  arc  operation   in  going  from  argon  to  nitrogen  to 

helium  plasmas.   Moreover,  the  increased  values  of  E    , 
c  x,a 

which  reflect  the  dissociation  energy  of  N„  and  the  rela- 
tively high  ionization  potential  of  He,  necessarily  result 
in  the  proportionately  higher  values  of  q   for  nitrogen 
and  helium  which  are  indicated  in  Figure  III. 9.   On  the 
basis  of  the  increased  values  of  total  wall  heat  flux  for 
nitrogen  and  helium,  it  might  be  expected  that  the  mean 
temperatures  of  these  two  gases  would  be  higher  than  the 
corresponding  values  of  T     for  argon.   However,  this  is 
clearly  not  the  case.   As  indicated  in  Figure  III. 13, 

the  T    values  for  nitrogen  and  helium  are  considerably 
m,a 

lower  than  the  T    values  for  argon.   In  fact,  at  I  =  50 

m,a 

amp,  the  experimentally  determined  value  of  T    for 
helium  is  more  than  an  order  of  magnitude  lower  than  the 
corresponding  value  for  argon,  and  yet  the  calorimetri- 
cally  measured  value  of  q   for  He  exceeds  the  correspond- 
ing  value  for  Ar  by  a  factor  of  nearly  six.   While  this 
behavior  is  seemingly  anomalous  to  the  rules  of  convection 
heat  transfer,  it  is  recognized  that  any  attempt  at  ex- 
planation on  the  basis  of  the  global  measurements  of  this 
study  would  be  fatuous.   While  a  partial  explanation  for 
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the  increased  wall  heat  transfer  for  the  helium  arc  over 
that  for  argon  can  be  inferred  from  the  characteristically 
higher  values  of  thermal  conductivity  for  He  (e.g.,  at 
T  =  5,000°K,  kRe  -  9  k^37'401),  a  complete  explanation 
cannot  be  gained  without  consideration  of  the  microscopic 
transport  phenomena  which  characterize  these  gases. 

The  order-of -magnitude  difference  between  the  mean 
temperatures  in  argon  and  helium  can  be  explained  on  the 
basis  of  the  relative  specific  heats  for  these  two  gases. 
That  is,  while  the  experimentally  determined  mean  en- 
thalpies for  He  and  Ar  are  nearly  equal  (Tables  IV. 1  and 
IV. 2),  the  specific  heat  for  helium  is  nearly  an  order  of 
magnitude  higher  than  the  corresponding  value  for  argon. 
When  evaluated  at  the  experimentally  determined  values 
of  h    for  argon  and  helium  (Tables  IV. 1  and  IV. 2)  ,  the 

III  f  cl 

ratio  of  the  constant  pressure  specific  heats  (He  to  Ar) 
is  9.1[38'391. 

IV. 4   Comparisons  for  the  Anode  Region 
Before  comparing  experimental  and  theoretical 
anode  heat  transfer  results,  it  is  necessary  to  discuss 
the  theoretical  model  and  its  relevance  to  measured 
quantities.   Since  the  major  portion  of  the  experimental 
data  acquired  in  this  study  is  for  argon,  the  following 
discussion  and  comparisons  will  deal  only  with  this  gas. 
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An  anode  heat  transfer  model  which  is  frequently 
found  in  the  literature  '  '  ~  is  one  based  upon  an 
energy  balance  over  the  surface  of  the  anode . 


<?r  +  c   <?e    <?abl   <?. 


**€ 


77777777777777777777 


JL 
T 


Figure  IV. 10   Anode  Heat  Balance  Model 


With  reference  to  Figure  IV. 10,  this  balance  is  of  the 
form 


•     •       •     • 


Q,  =  Q    +Q  -  Q  i»i  -  Q   / 
WA   wr+c    e   *abl    e  ' 


(IV. 18) 


where 


Q    =  rate  of  energy  loss  to  the  anode  cooling 
water , 
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Q  +  =  rate  of  energy  transfer  to  the  anode  sur- 
face by  convection  and  radiation  from  the 
plasma  column, 

Q  =  rate  of  energy  transfer  to  the  anode  sur- 
face by  electrons, 

Q  ,  ,  =  rate  of  energy  transfer  from  the  anode  sur- 
face due  to  ablation  of  the  anode  material, 
and 


Q    =  radiation  emissive  power  of  the  anode  sur- 
face. 


On  the  basis  of  observed  negligible  degradation  of 
the  anode  surface  after  some  60  hours  of  operation  and 
from  calculations  for  a  black  surface  emitting  at  the 
melting  temperature  of  copper,  the  last  two  terms  in  Equa- 
tion IV. 18  are  known  to  be  insignificant  compared  to 
calorimetrically  measured  values  of  Q,  and  can  therefore 
be  neglected.   In  addition,  the  energy  transferred  to 
the  anode  surface  by  the  electrons  is  commonly  obtained 

r  6 1 

from  the   expression 

Qe  =  IU    <HkTe  +  V    +  *A!'  <IV-19) 
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where 

I  5 

—  j  ^T   =  thermal  energy  of  the  electrons  as  they 

enter  the  region  of  "anode  fall", 

IU-      =  additional  energy  acquired  by  the  elec- 
trons due  to  acceleration  through  the 
anode  fall  region, 

I<|>_     =  "condensation"  energy  of  the  electrons 
penetrating  the  anode  surface  (<J>A  = 
work  function  of  the  anode  material) ,  and 

£  =  fraction  of  energy  retained  by  the  elec- 
trons after  collisions  with  heavy  parti- 
cles within  the  anode  fall  region. 

In  addition  to  assuming  negligible  Qahl  and  Q  , 

f  6  71 
Pfender,  et  aL   '   ,  in  working  with  Equations  IV. 18  and 

IV. 19  assume  a  value  of  unity  for  £  and  assume  that  U_ 
may  be  represented  by  the  difference  between  the  anode 
potential  (measured  with  respect  to  the  cathode)  and  the 
potential  obtained  by  linearly  extrapolating  the  fully- 
developed  potential  gradient  (E    )  to  the  location  of 

x  /  a 

the  anode  (see  Figure  IV. 11).  In  using  these  assumptions 
with  selected  values  of  T  =  10,000°K  and  <j>  =  3.7  volts, 
Equations  IV. 18  and  IV. 19  reduce  to 
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QA  =  Qr+c  +  X    (5'85  +  V*  (IV. 20) 

If  the  anode  is  preceded  by  a  thermally  fully-developed 

region  and  if  the  mean  enthalpy  of  the  plasma  at  the  anode 

remains  close  to  h    ,  Pfender,  et  al. argue  that  Q 

m,a  3         r+c 

becomes  the  asymptotic  value  of  the  total  wall  heat  flux 
(q  ) .   Accordingly,  Equation  IV. 20  takes  the  form 

CL 

QA  =  qa  +  26.6  I  (5.85  +  UA) ,  (IV. 21) 

where  the  constant  26.6  accounts  for  the  area  of  the  anode 

2 

surface  (Q   and  qa  have  units  of  [Btu/f t  -min] )  and  the 

Btu 

appropriate   conversion    factor    (0.0569   — rr ■ — ). 

watt-mm 

Using  the  experimental  data  of  this  study  for  q 

a 

and  \J    ,  a  comparison  is  made  between  the  values  of  QA 
calculated  from  Equation  IV. 21  and  the  measured  values  of 
this  study  (Figure  III. 15).   The  results  of  this  comparison 
are  shown  in  Table  IV. 3,  wherein  it  is  noted  that  the 
theoretical  model  overpredicts  the  corresponding  experi- 
mental values  by  a  factor  of  approximately  two.   It  is  felt 
that  a  portion  of  this  rather  gross  overprediction  is  due 
to  errors  associated  with  assuming  that  K   =  1.0  and 
assuming  that  U,  can  be  determined  by  the  aforementioned 
linear  extrapolation  procedure.   These  errors  are  made 
evident  by  the  anode  energy  balance  measurements  of  this 
study  which  reveal  a  significant  increase  (~  30%  above 
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the  value  of  h    )  in  the  enthalpy  of  the  gas  as  it  passes 
through  the  anode  region.   In  order  to  obtain  such  an 
increase  in  enthalpy,  there  must  occur  a  significant  trans- 
fer of  electrical  energy  to  the  gas  in  excess  of  that 

required  to  maintain  the  asymptotic  conditions  (h    and 

-1  m,a 

thus  q  ) .   The  electrons,  therefore,  will  in  fact  lose  a 
a 

portion  of  their  energy  (IU_)  to  the  heavy  particles  in 
the  fall  region  and,  consequently,  E,   will  be  considerably 
less  than  unity.   The  failure  of  the  theoretical  model 
(Equation  IV. 21),  when  applied  to  the  anode  configuration 
of  the  constricted  tube  arc,  to  account  for  this  additional 
transfer  of  electrical  energy  is  perhaps  best  understood 
by  considering  a  typical  anode  region  potential  distri- 
bution.  With  reference  to  Figure  IV. 11,  the  theoretical 
model  assumes  that,  of  the  total  anode  potential  rise,  E  , 
the  portion,  E  -U  ,  is  required  to  maintain  the  asymptotic 

value  of  h    and  the  remaining  portion,  U  ,  results  in 
m ,  a  a 

energy  transfer,  via  the  electrons,  to  the  anode  surface 
(i.e.,  IU  ) .   In  contrast,  the  experimental  results 
indicate  that,  in  addition  to  that  portion  of  the  total 
anode  potential  rise  required  to  maintain  h    ,  a  portion 
of  E,  is  also  needed  to  sustain  the  increase  in  the 
enthalpy  of  the  gas,  Ah_ ,  as  it  passes  through  the  anode 
region.   Consequently,  a  fraction  of  the  energy  acquired 


Ill 


Figure  IV. 11   Partitioning  of  Total  Anode  Fall 
Potential 


by  the  electrons  (IUA)  is  in  fact  transferred  to  the  gas 
through  electron-heavy  particle  collisions  and  cannot 
contribute  directly  to  the  total  anode  heat  transfer. 
Hence,  it  is  concluded  that  the  overpredictions  for  Q 
which  result  from  the  anode  heat  balance  model  are  at 
least  partially  due  to  the  combination  of  assuming  that  £ 
is  unity  and  that  the  anode  fall  potential  is  given  by 
the  linear  extrapolation  procedure. 

To  further  consider  these  assumptions  as  the  major 
sources  of  error  in  Equation  IV. 21,  the  data  of  this  study 
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were  used  to  determine  both  the  potential  drop  required 
to  sustain  the  gas  enthalpy  rise  Ah  and  the  potential 
drop,  to  be  designated  U  ,  ,  which  results  in  the  transfer 
of  electron  energy  directly  to  the  anode  surface  (ID  ,  ) . 
In  this  connection  and  with  reference  to  Figure  IV. 11,  it 
is  assumed  that  the  total  anode  potential  rise,  E  ,  can  be 
subdivided  into  the  fraction  of  this  potential  required 
to  sustain  the  asymptotic  conditions,  E-U./E. ;  the  frac- 

r\    r\    A. 

tion  of  this  potential  required  to  sustain  the  anode  re- 
gion enthalpy  rise  (Ah  ) ,  U  -UA/EA  ;  and  the  fraction  of 
this  potential  resulting  in  energy  transfer  to  the  anode 
surface  by  the  electrons,  Ua,  /E  .   According  to  this  sub- 
division,  and  retaining  the  values  of  T   and  4>   used  in 
Equation  IV. 21,  energy  balance  considerations  require  that 


Q   =  q   +  26.6  I  (5.85  +  Unl  )  ,         (IV. 22) 
A    a  a 


and  that 

I  (UA  -  UA,  )  =  rfiAhA  ,  (IV.  23) 

where,  again,  Ah   is  the  enthalpy  acquired  by  the  gas  in 
traversing  the  anode  region.   The  anode  heat  flux  predic- 
tions which  result  from  calculations  performed  with  this 
modified  theoretical  model  are  presented  in  Table  IV. 4 
along  with  the  experimentally  determined  values  of  AhA« 
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While  the  values  of  U  ,  computed  from  Equation  IV. 23 

are  significantly  lower  than  those  given  by  the  linear 

extrapolation  method  (U  ) ,  their  use  in  Equation  IV. 22, 

along  with  experimental  values  of  q  ,  still  result  in 

a 

sizeable  overpredictions  for  the  total  anode  heat  flux. 
It  appears  then  that  more  appropriate  values  for  T  and 
<J>  must  be  determined  in  order  to  bring  the  theoretical 
predictions  into  better  agreement  with  experiment.   These 
additional  considerations  would  necessarily  require  further 
knowledge  of  electron  temperature  distributions  in  the 
vicinity  of  the  anode  as  well  as  a  more  rigorous  consid- 
eration of  the  plasma-electrode  surface  interaction.   Ac- 
cordingly, the  use  of  this  analytical  model  to  predict  the 
total  heat  transfer  to  the  anode  of  a  constricted  tube 
arc  appears  to  be  of  questionable  value.   This  is  espec- 
ially true  when  the  model  is  compared  with  existing  em- 
pirical correlations,  such  as  those  appearing  in  Figure 
III. 15.   In  connection  with  these  empirical  correlations, 
it  is  noted  that  in  cases  where  the  anode  is  not  preceded 
by  a  fully-developed  region,  it  is  necessary  to  include 

a  dependency  on  arc  length  and  operating  gas  flow  rate, 

[29] 
as  is  the  case,  for  example,  of  the  Eberhart  and  Seban 

anode  heat  transfer  correlation. 
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CHAPTER  V   SUMMARY  AND  CONCLUSIONS 

This  work  provides  experimental  results  for  the 
axial  distribution  of  the  electric  field  intensity  and 
total  wall  heat  flux  for  the  inlet,  asymptotic,  and  anode 
regions  of  atmospheric  argon,  helium,  and  nitrogen  arc 
plasmas.   Whereas  data  are  presented  for  a  wide  range  of 
laminar  flow  rates  (0.03-0.11  lb/min)  and  arc  currents 
(35-240  amp)  in  argon,  the  amount  of  data  presented  for 
helium  and  nitrogen  arcs  is  considerably  less  extensive. 
The  range  of  operation  for  these  two  gases  was  restricted 
by  the  power  supply  employed  in  this  study.   In  addition 
to  the  development  of  the  diagnostics  required  for  the 
electric  field  intensity  and  total  wall  heat  flux  measure- 
ments, considerable  effort  was  given  to  the  development 
of  a  unique  gage  to  measure  the  radiative  component  of 
the  total  wall  heat  flux.   This  instrument  has  been  used 
successfully  in  providing  the  radiative  wall  heat  flux 
data  which  are  presented  for  the  asymptotic  region  in 
argon.   The  accuracy  of  the  experimental  measurements, 
as  well  as  the  consistency  of  the  arc  facility,  is  assessed 
on  the  basis  of  the  observed  precision  of  local  energy 
balance  checks.   Further  confidence  in  the  experimental 
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measurements  is  inferred  from  the  favorable  agreement 
which  exists  between  the  data  of  this  study  and  the  experi- 
mental results  reported  by  other  investigators. 

Through  comparisons  between  the  experimental  data 
and  available  equilibrium  and  nonequilibrium  numerical 
solutions,  an  attempt  is  made  to  ascertain  the  adequacy  of 
single-fluid  models  to  accurately  describe  overall  arc 
energetics.   On  the  basis  of  the  measurements  themselves 
and  the  comparison  with  theoretical  results,  the  following 
conclusions  were  derived. 

1.  The  increase  in  duct  wall  potential  with  current 
drawn  from  the  segment  is  primarily  a  "cold" 
boundary  layer  effect  and,  therefore,  the 
"floating"  duct  wall  potential  provides  a  more 
accurate  measure  of  the  actual  plasma  potential, 
at  least  for  the  copper  wall-argon  plasma 
combination  employed  in  this  study. 

2.  Variations  in  the  duct  entrance  conditions  can 
have  a  significant  effect  on  the  experimentally 
determined  and  numerically  predicted  electric 
field  intensity  and  total  wall  heat  flux  dis- 
tributions for  the  inlet  region  of  the  con- 
strictor.  As  a  result,  meaningful  comparisons 
between  experiment  and  theory  can  be  made  only 
for  the  asymptotic  region. 
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3.  Closed  form  analytical  solutions  which  are 
based  upon  the  use  of  linearized  property 
relationships  and  which  assume  negligible 
radiation  are  inadequate  for  predicting  electric 
field  intensity  and  total  wall  heat  flux.   This 
is  particularly  evident  at  the  higher  current 
levels  in  argon  and  nitrogen. 

4.  To  within  the  combined  experimental  and  theoreti- 
cal uncertainties,  the  numerical  solution  of  a 
rigorous  single-fluid  model  is  suitable  for  use 
in  predicting  the  electric  field  intensity  and 
total  wall  heat  flux  for  the  asymptotic  region 

in  argon,  at  least  over  the  range  of  conditions 
investigated  in  this  study.   The  equilibrium 
solution  does  not,  however,  appear  suitable 
for  predicting  the  asymptotic  mean  temperature 
or  the  radiative  wall  heat  flux  over  the  same 
range  of  conditions  in  argon.   There  is  a  small, 
but  discernable,  nonequilibrium  effect  at  the 
lower  arc  currents  in  argon,  as  a  result  of 
which  the  equilibrium  solution  tends  to  over- 
predict  both  the  electric  field  intensity  and 
total  wall  heat  flux  for  currents  less  than 
100  amp.   Moreover,  the  tendency  of  the  equili- 
brium solution  to  underpredict  these  two  para- 
meters at  arc  currents  in  argon  greater  than 
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100  amp  is  believed  to  be  due  to  a  chemical 
nonequilibrium  effect  which  renders  the 
equilibrium  solution  deficient  in  providing 
for  sufficient  energy  loss  to  the  tube  walls 
by  conduction. 

5.  Solutions  which  properly  account  for  thermo- 
chemical  nonequilibrium,  as  well  as  for  all 
pertinent  transport  phenomena,  provide  pre- 
dictions for  the  electric  field  intensity  and 
the  convective  and  radiative  wall  heat  transfer 
which  are  in  excellent  agreement  with  the 
experimental  results  over  the  entire  range  of 
asymptotic  conditions  investigated  in  argon. 

6.  Atmospheric  helium  arcs  are  characterized  by  a 
significant  departure  from  local  thermodynamic 
equilibrium  (at  least  for  arc  currents  less 
than  250  amp)  and,  as  a  result  of  these 
departures,  single-fluid  models  are  totally 
inadequate  in  representing  their  thermal  state. 

With  reference  to  the  last  three  of  the  above  con- 
clusions, it  is  presumed  that  further  comparisons  between 
experiment  and  equilibrium  and  nonequilibrium  solutions, 
if  extended  to  include  arc  conditions  characterized  by 
higher  order  nonequilibrium  effects,  would  lead  to  an 
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improved  understanding  of  nonequilibrium  arc  physics. 
Accordingly,  it  is  recommended  that  the  present  arc 
facility  be  extended  to  permit  stable,  low  current 
(I  <  35  amp)  operation  in  argon  and  to  permit  stable 
operation  over  a  greater  range  of  operating  conditions 
in  helium  and  nitrogen.   It  is  also  recommended  that  the 
solution  by  Bower  and  Incropera     be  expanded  to  provide 

equilibrium  predictions  for  nitrogen  arcs  and  that  the 

[251 
solution  by  Clark     be  appropriately  extended  to  afford 

nonequilibrium  arc  predictions  for  helium  and  nitrogen. 

In  addition  to  the  above  conclusions  for  the  inlet 
and  asymptotic  regions,  further  comparisons  between  exper- 
iment and  theory  for  the  anode  region  indicate  that  the 
large  (50-100%)  overpredictions  for  the  total  anode  heat 
flux  which  result  from  use  of  a  heat  balance  model  are 
primarily  due  to  the  relatively  large  uncertainties  in- 
volved in  specifying  associated  values  of  the  physical 
parameters  (U  ,  T  ,  and  <j>  )  .   Consequently,  the  theoreti- 
cal  model  is  thought  to  be  ill-suited  for  providing  de- 
sign estimates  for  the  total  heat  flux  to  the  anode  of 
the  arc  constrictor,  particularly  in  comparison  with 
estimates  which  are  based  upon  empirical  correlations.   In 
this  respect,  further  experimental  and  theoretical  study 
of  plasma-electrode  surface  interactions  is  recommended. 
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APPENDIX   RADIATION  GAGE 

A.l   General  Remarks  and  Review  of  Experimental 
Studies  Relating  to  Radiative  Heat  Flux  Measurements 

Since  thermal  radiation  is  often  a  primary  mode 
of  energy  transfer  in  a  constricted  arc  plasma,  knowledge 
of  the  radiative  contribution  to  the  total  wall  heat  flux 
is  of  concern  in  both  the  engineering  application  and 
the  laboratory  use  of  plasma  arc  generators.   In  addition, 
knowledge  of  the  radiative  wall  heat  flux  is  of  value  in 
obtaining  an  improved  understanding  of  arc  physics  and 
plasma  flow  modeling.   Accordingly,  a  portion  of  this 
study  was  directed  towards  the  design,  fabrication,  and 
use  of  appropriate  instrumentation  to  measure  the  radia- 
tive component  of  the  total  heat  flux  to  the  wall  of  an 
arc  constrictor. 

While  there  are  numerous  studies  in  the  literature 
dealing  with  the  measurement  of  plasma  radiative  proper- 
ties ,  the  author  knows  of  no  experimental  attempt  to 
deal  explicitly  with  the  measurement  of  the  radiative 
and  convective  components  of  the  total  wall  heat  flux  for 
an  arc  constrictor.   There  are  described  in  the  litera- 
ture, however,  two  diagnostic  techniques  which  would  appear 
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applicable  to  this  particular  measurement.   The  first  of 
these  involves  a  total  radiation  detector  which  utilizes 
a  thin-film  resistance  sensor.   This  sensor  consists  of 
a  thin  metallic  strip  (usually  platinum)  which  is  depo- 
sited on  a  substrate  material  (pyrex,  quartz,  etc.)  and 
has  been  used  primarily  for  transient,  short  duration 
measurements.   The  metallic  strip  acts  as  a  passive  re- 
sistance element  in  an  electrical  circuit  designed  to 
record  the  variation  of  current  through  the  sensor  upon 
its  exposure  to  a  step  input  of  radiation.   From  know- 
ledge of  the  temperature  coefficient  of  resistance  of  the 
metallic  element  and  the  thermal  conductivity  of  the  sub- 
strate, the  heat  flux  to  the  substrate  can  be  inferred 
by  analyzing  the  time  variation  in  the  surface  temperature 
as  indicated  by  the  electrical  output  from  the  sensor. 
Further  discussion  of  this  diagnostic  technique  is  provided 

r  6 1 1 

by  Reece  and  Theoclitus     .   A  major  advantage  of  this 

method  is  that  the  calibration  and  data  reduction  procedures 

[62,63] 
are  relatively  straight  forward  and  well  established 

While  the  time  limit  of  sensitivity  (saturation 

time)  for  the  thin-film  sensor  is  typically  on  the  order 

of  100  msec,  this  technique  can  be  modified  for  use  with 

steady  state  radiation  sources.   As  discussed  by  Hardy 

and  Paddock  *-    ,  this  feature  is  made  possible  by  placing 
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a  fast-acting  shutter  mechanism  between  the  radiation 
source  and  the  sensor.   This  modification  would,  of  course, 
be  required  if  the  thin-film  detector  were  to  be  used  for 
the  radiation  measurements  of  this  study.   However,  here- 
in lies  the  major  disadvantage  of  this  method.   In  parti- 
cular, the  imposition  of  interface  components  (shutter 
mechanism,  etc.)  would  place  the  sensor  sufficiently  far 
from  the  plasma  source  to  require  detailed  consideration 
of  the  spectral  sensitivity  of  the  resulting  cavity  walls. 

A  second  device  has  been  employed  by  researchers 
at  Harvard  to  measure  the  total  radiation  from  a  steady 
state,  constricted  tube  arc*   This  device  uses  a  re- 
cessed sensing  element,  across  which  a  temperature 
difference  is  induced  which  is  proportional  to  the  ab- 
sorbed radiation.   The  sensing  element  consists  of  a 
blackened  constantan  shim  which  is  positioned  in  a  cir- 
cumferential (with  respect  to  the  axis  of  the  cylindrical 
duct)  slot.   The  slot  walls  are  silver  coated  in  order  to 
reflect  a  greater  percentage  of  radiation  to  the  sensor, 
and  the  sensing  element  is  recessed  1.5  inches  in  an 
attempt  to  minimize  convective  contributions  to  the 
measured  heat  flux.   Calibration  of  the  radiation  meter 
was  attempted  by  comparing  the  detector  output  against 


*  A  comparison  is  made  in  Chapter  III  between  the  Harvard 
results  and  the  radiation  data  obtained  in  this  study. 
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several  radiation  standards,  which  include  filament  and 
gas-filled  lamps  and  globar  units. 

This  diagnostic  technique  is  subject  to  several 
sources  of  error  for  which  compensation  cannot  be  entirely 
made  by  the  calibration  procedure.   Consequently,  accurate 
use  of  this  device  is  contingent  upon  the  successful 
isolation  and  evaluation  of  these  sources  of  error  and  the 
subsequent  introduction  of  appropriate  correction  factors. 
This  matter  deserves  further  discussion,  for  it  serves  to 
illustrate  the  problems  associated  with  radiative  heat 
flux  measurements,  particularly  with  techniques  which 

employ  recessed  sensors. 

[12] 
As  discussed  by  Emmons     ,  certain  correction 

factors  must  be  applied  to  the  experimental  results  ob- 
tained by  use  of  the  aforementioned  device.   This  is 
necessary  in  order  to  account  for  the  geometrical  differ- 
ences between  the  conditions  of  calibration  and  actual 
measurement.   Specifically,  account  must  be  taken  of  the 
difference  between  the  relative  positions  (with  respect 
to  the  sensor  location)  of  the  radiation  standard  and  the 
arc.   In  addition  the  variation  in  geometry  which  results 
from  the  fact  that  the  calibration  standard  (globar  unit) 
is  a  surface  emitter,  whereas  the  arc  is  a  volume  emitter, 
must  be  considered.   These  correction  factors  may  be 
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obtained  from  geometrical  considerations  relating  to  the 
calibration  and  actual  measurement  configurations .   How- 
ever, aside  from  the  geometrical  effects,  there  remains 
a  possible  source  of  error  due  to  the  spectral  sensitivity 
of  the  component  parts,  particularly  the  slot  walls.   If 
the  reflectivity  of  the  slot  walls  varies  appreciably  over 
the  wavelength  range  of  the  incident  plasma  radiation,  and 
if  this  variation  is  to  be  incorporated  within  the  cali- 
bration of  the  detector,  then  the  radiation  standard  must, 
of  course,  have  the  same  spectral  intensity  distribution 
as  the  plasma.   Unfortunately,  such  calibrated  sources 
are  not  available. 

Another  deficiency  inherent  in  this  radiation  gage 
is  the  fact  that  the  higher  plasma  radiation  intensities 

could  only  be  determined  from  a  linear  extropolation  of 

r  4  7 1 

the  calibration  data     .   That  is,  the  detector  output 

during  actual  measurement  was  a  factor  of  approximately 
six  higher  than  the  maximum  calibrated  detector  output 
(as  limited  by  the  burnout  temperature,  1625°C,  of  the 
globar  unit) .   Techniques  employing  recessed  sensors  are 
therefore  calibration  limited  in  the  sense  that  if  the 
possible  variations  in  the  spectral  response  of  the  cavity 
walls  (as  well  as  the  full  output  range  of  the  detector) 
are  to  be  incorporated  within  the  detector  calibration, 
the  calibration  source  must  have  the  same  radiation 
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intensity,  in  both  magnitude  and  spectral  distribution, 
as  the  plasma  for  which  measurements  are  desired. 

A. 2   Radiative  Heat  Transfer  Gage  -  General 
Design  and  Operation  Considerations 

From  evaluation  of  the  aforementioned  measurement 
techniques,  three  primary  constraints  were  initially 
imposed  upon  the  proposed  detector  design.   First,  the 
detecting  unit  must  not  alter  the  geometry  of  the  con- 
strictor and,  in  particular,  the  use  of  recessed  sensing 
elements  should  be  avoided.   Secondly,  the  detector  must 
not  be  calibration  limited  in  the  sense  that  generation 
of  its  output  characteristic  is  dependent  upon  an  initial 
calibration  against  a  standard  radiation  source.   Thirdly, 
the  presence  of  the  radiation  gage  should  not  significantly 
alter  local  conditions  in  the  plasma.   This  is  necessary 
if  useful  comparisons  between  theory  and  measurements 
are  to  be  realized. 

To  adhere  to  the  above  criteria ,  it  seemed  rea- 
sonable to  attempt  to  alter  the  wall  surface  characteris- 
tics to  permit  isolation  and  measurement  of  the  radiative 
heat  flux.   The  use  of  a  transparent  "window"  affords  a 
convenient  means  of  isolation;  however  presently  avail- 
able optical  materials  cannot  provide  uniform  trans- 
mission over  the  entire  spectral  range  of  possible  plasma 
radiation.   In  addition,  the  window  would  have  to  be 
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thermally  isolated  from  the  sensor  so  as  not  to  contri- 
bute to  the  measured  heat  flux  through  conduction. 

Further  insight  into  a  possible  technique  for  de- 
lineating between  the  radiative  and  convective  components 
to  the  total  wall  heat  flux  was  gained  from  consideration 
of  the  optical  characteristics  of  the  duct  wall  surface. 
Expressing  the  radiative  component  in  terms  of  the  mono- 
chromatic absorptance  (a,)  and  the  plasma  irradiation 


(H,  ) ,  the  total  wall  heat  flux  takes  the  form 


^conv       ^rad 


=  q  +   /   cuH,dX  (A.l) 

^conv        ;       X    X 
o 


4 


The  emissive  power  of   the  duct  wall  surface    (ecrT      )    is 

W 

neglected,  since  its  magnitude  is  small  in  comparison  to 

the  other  terms  of  the  equation.   Assuming,  for  the  moment, 

that  a,  is  known  and  that  q  may  be  determined  by  calori- 

metric  measurement,  there  remain  two  unknowns  in  Equation 

A.l,  q     and  H, .   The  appearance  of  these  two  unknowns 
'  ^conv      X        rtr 

suggests  the  use  of  two  surfaces  (which  differ  in  spectral 
sensitivity) ,  in  which  case  the  application  of  Equation 
A.l  for  surface-1  yields 
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00 

I  ql  =  qoonv,l  +  /aX,l  HX,1  dX  (A-2) 


and,  for  surf ace-2 , 


q2  =  ^conv,2  +  la\,2   HX  ,2  dX  (A'3) 


If  both  surfaces  are  located  at  the  same  axial  location 
in  the  duct,  it  may  be  assumed  that  both  the  incident 
irradiation  and  the  convective  heat  transfer  to  each  sur- 
face are  the  same  (i.e.,  H,  ,  =  H.    =  H,  and  q     .  = 

'      A,l     \  ,2  X  ^conv,l 

qco   2  =  qconv^ *   Under  these  assumptions,  Equations 

A. 2  and  A. 3  can  be  equated  on  the  basis  of  equal  convection 

to  give 


00  00 

ql  -  q2  =  K,!  HX  dX  -  /aX,2  HX  dX  •        (A-4) 

0    '  0    '  / 


To  achieve  the  largest  possible  difference  between 
the  calorimetrically  measured  values  of  the  total  wall 
heat  flux  (q.  -  q2)  and  thus  to  realize  the  maximum 
experimental  resolution,  the  difference  between  the  radia- 
tion absorption  terms  in  Equation  A. 4  should  be  maximized. 
This  suggests  the  use  of  a  highly  absorbing  (black)  sur- 
face in  combination  with  a  highly  reflecting  surface. 
For  the  optically  black  surface,  say  surf ace-1 ,  it  is 
assumed  that  suitable  coating  materials  are  available  for 
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which  the  monochromatic  absorptance  is  independent  of 
wavelength  over  the  entire  spectral  range  of  the  plasma 
irradiation.   For  such  a  surface,  the  radiative  wall 
heat  flux  component  reduces  to 

oo  00 

!  a\   I  H\  dX  =   ax  /Hx  dX  =  aiH  •  (A-5> 

0     '  0 

Similarly  for  the  reflecting  surface  (surface-2) , 

OO  CO 

!    ax    2      EX    dX    =   a2    ^  dX  =  a2H  *  (A,6) 

0     '  0 

With  these  additional  assumptions,  Equation  A. 4  reduces 
to 

q-j^  -  q2  =  H(a1  -  a^)  (A. 7) 

Knowledge  of  the  absorptance  values,  a,  and  a_,  along  with 
the  experimentally  determined  values  of  q,  and  q2  will 
permit  evaluation  of  the  plasma  irradiation.   This  value 
of  H  can  then  be  employed  in  either  of  the  following 
expressions  to  yield  the  value  of  the  convective  wall  heat 
transfer : 

q,  =  q      +  a  H  (A. 8) 

1    conv     1 

or, 

q0  =  q      +  a0H  .  (A. 9) 

^2    ^conv    2 
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To  determine  the  radiative  wall  heat  flux  to  the  standard 

(copper)  duct  segments,  this  value  of  q     is  subtracted 

conv 

from  the  total  wall  heat  flux  results  (obtained  at  the 
same  operating  conditions  and  axial  location)  determined 
from  previous  calorimetric  measurements . 

As  shown  in  Figure  A.l,  the  radiation  gage  consists 
of  two  semicylindrical  copper  segments  which  are  thermally 
isolated  from  one  another,  individually  water  cooled, 
and  supplied  with  individual  AT-transducers .   An  optical 
black  (Nextel  101-C10  by  Minnesota  Mining  and  Manufactur- 
ing Company)  coating  is  applied  (and  then  baked)  to  the 
interior  cylindrical  surface  of  one  of  the  half-segments, 
while  the  opposing  half-segment  is  coated  with  a  highly 
reflecting,  thermally  evaporated  aluminum  film.*   The 
overall  dimensions  of  the  assembled  radiation  gage  are 
precisely  the  same  as  those  of  a  standard  duct  segment. 

A. 3  Discussion  of  the  Parameters  and 
Assumptions  Affecting  the  Operation 
of  the  Radiation  Gage 
While  the  above  radiation  gage  design  appears  to 
satisfy  all  of  the  initially  posed  constraints  and 
specifications ,  the  following  two  basic  assumptions  re- 
quire further  discussion: 


*  The  fabrication  and  optical  properties  of  the  reflecting 
surface  are  discussed  in  Section  A. 4. 
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1.  The  presence  of  the  radiation  gage  does  not 
significantly  alter  the  thermal  state  of 
the  plasma  in  the  vicinity  of  the  gage  and, 

2.  The  reflecting  surface  is  characterized  by  a 
constant  value  of  absorptance,  a?,    over  the 
entire  spectral  range  of  the  incident  plasma 
irradiation  (H, ) . 

The  validity  of  the  first  assumption  is  critical 
to  the  success  of  this  diagnostic  technique.   That  is, 
if  the  pressence  of  the  highly  absorbing  and  reflecting 
surfaces  results  in  a  significant  change  in  the  local 
thermal  state  of  the  plasma,  then  the  assumptions  of 
equal  convection  and  irradiation  required  by  Equation 
A. 7  would  be  rendered  questionable.   Moreover,  any  such 
change  in  the  local  thermal  plasma  state  would  preclude 
combining  the  convective  wall  heat  flux  (determined  from 
use  of  the  radiation  gage)  with  the  total  wall  heat  flux 
data  (obtained  from  calorimetric  measurements  employing 
the  standard,  copper-wall  duct  segments)  to  determine 
the  radiation  flux. 

Wall  heat  transfer  measurements  were  made  at  the 
radiation  gage  segment  and  the  eight  adjoining  standard 
segments  (four  upstream  and  four  downstream)  in  an  attempt 
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to  experimentally  determine  if  there  existed  any  effect 
on  the  plasma  radiosity  distribution  that  could  be 
attributed  to  the  spectral  characteristics  of  the  gage 
surfaces.   This  experimental  investigation  in  the 
asymptotic  region  of  an  argon  arc  was  performed  using  both 
black-versus-reflector  and  ref lector-versus-ref lector 
surface  combinations .   The  results  showed  no  discernable 
difference  between  the  total  wall  heat  flux  to  the  seg- 
ments adjoining  the  radiation  gage  and  the  q  values  which 

a 

were  measured  in  the  absence  of  this  gage.   Moreover,  the 
total  wall  heat  flux  to  both  of  the  reflecting  surfaces 
for  the  case  of  the  ref lector-versus-ref lector  combina- 
tion compares  favorably  with  the  measured  total  wall 
heat  flux  to  the  reflecting  surface  of  the  black-versus- 
reflector  combination.*   During  the  course  of  these 
experiments ,  the  duct  wall  potential  distribution  was 

monitored  and,  as  in  the  case  of  q  ,  there  were  no  dis- 

a 

cernable  variations  in  this  parameter.   Accordingly,  it  is 

concluded  that,  to  within  the  experimental  uncertainty  of 

q  (±6%)  and  E    (±  3%)  measurements,  the  presence  of  the 
a  x ,  a 

radiation  gage  does  not  significantly  alter  the  local 
thermal  state  of  the  plasma. 


*  These  observations  are  discussed  further  in  Section  A. 5. 
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With  respect  to  the  second  of  the  two  aforementioned 
assumptions,  consideration  must  be  given  to  the  spectral 
distribution  of  the  incident  plasma  irradiation,  H  =  H(A), 
before  a  suitable  reflecting  material,  for  which  a2  ^  a2(X) 
over  the  spectral  range  of  H, ,  can  be  selected.   The  fol- 
lowing discussion  will  be  restricted  to  the  case  of  an 
argon  plasma  since  the  radiative  wall  heat  flux  data  ob- 
tained in  this  study  pertain  only  to  this  gas. 

r  66 1 

On  the  basis  of  the  theoretical  work  of  Yakubov     , 
the  emission  spectrum  from  an  atmospheric  argon  plasma 
is  known  to  cover  a  wavelength  band  extending  from  the 
vacuum  ultraviolet  to  the  infrared.   Of  particular  interest 

o 

to  this  discussion  is  the  vacuum  ultraviolet  (A  <  2000  A) 
portion  of  the  emission  spectrum  which  results  from  bound- 
and  free-bound  transitions  to  the  ground  state  (resonance 
lines  and  ionization  continuum,  respectively) .   As  indi- 
cated by  Yakubov,  and  on  the  basis  of  the  optical  depth 

and  radiative  recombination  rate  calculations  performed 

[671 
by  Horn     ,  the  radiation  emitted  from  these  transitions 

is  strongly  reabsorbed,  and,  therefore,  does  not  contri- 
bute significantly  to  the  balance  of  energy  at  the 
boundary  of  the  plasma.   Accordingly,  it  is  assumed  that 
there  are  no  significant  contributions  to  the  total 
plasma  radiation  incident  on  the  tube  wall  at  wavelengths 
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less  than  2000  A.   This  assumption  is  supported  by  the 
theoretical  calculations  of  Maecker  "-   ^  and  is  consistent 
with  the  plasma  radiation  model  employed  in  the  non- 

r  251 

equilibrium  solution  by  Clark    J  as  well  as  with  the 
experimental  radiation  source  density  employed  in  the 
equilibrium  solution  by  Bower  and  Incropera  *■    .   It 
remains,  however,  to  select  a  reflecting  surface  for  which 
the  monochromatic  absorptance  is  invariant  for  wavelengths 

o 

greater  than  2000  A. 

A. 4   Selection,  Fabrication,  and  Optical 
Characteristics  of  the  Reflecting  Surface 
While  the  spectral  sensitivity  of  a  metallic  sur- 
face is  enhanced  through  polishing,  the  optical  charac- 
teristics of  polished  metal  surfaces  are  both  difficult 
to  reproduce  and  to  maintain.   A  possible  alternative 
is  to  adjust  the  optical  characteristics  of  the  surface 
by  application  of  a  thin  metallic  coating  or  thin  film. 
Commonly  used  methods  of  thin  film  preparation  include 
thermal  evaporation,  sputtering,  and  electroplating,  or 
chemical  deposition.   The  electroplating  process  produces 
films  which  are  relatively  porous,  rough,  and  impure. 
This  process  is  therefore  generally  ill-suited  for 
applications  requiring  reproducible,  optical  quality  sur- 
faces.  Similarly,  film  deposition  by  sputtering  is  a 
relatively  slow  process  resulting  in  impure  and  low 
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density  coatings  which  are  unsuitable  for  use  as  ultra- 
violet reflectors.   Techniques  employing  thermal  eva- 
poration in  a  vacuum,  however,  have  proven  to  be  most 
satisfactory  for  precise  deposition  of  reproducible, 
optical  quality,  thin  film  coatings.   The  primary  con- 
trol variables  for  this  process  are  the  rate  of  deposi- 
tion and  the  evaporator  chamber  pressure. 

While  many  metal  films  (e.g.,  aluminum,  silver, 
gold,  and  rhodium)  are  highly  reflecting  at  wavelengths 
in  the  visible  and  infrared,  only  aluminum  films  main- 
tain a  high  reflectance  into  the  ultraviolet.   Through 
proper  preparation  and  aging,  aluminum  films  will  re- 
main better  than  90%  reflecting  at  wavelengths  greater 

o 

than  2000  A.   For  this  reason,  the  aluminum  thin  film  is 
particularly  suited  for  the  reflecting  surface  of  the 
radiation  gage. 

To  achieve  and  to  maintain  the  highest  possible 
reflectance  for  the  aluminum  film,  precise  control  is 
required  over  the  following  variables. 

1.  Film  purity. 

2.  Rate  of  deposition. 

3.  Evaporator  chamber  pressure. 

4.  Substrate  temperature. 

5.  Film  thickness. 

6.  Aging  conditions. 
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The  effect  of  each  of  these  variables  on  the 
reflectance  characteristics  of  evaporated  aluminum  films 

r  g  g  "I 

has  been  investigated  by  Madden     .   The  following 
discussion  provides  a  summary  of  this  work. 

The  purity  of  the  aluminum  has  a  slight  effect  on 
the  optical  characteristics  of  the  film  at  ultraviolet 
wavelengths,  with  purer  films  resulting  in  a  higher  percent 
reflectance.   There  is  little  improvement,  if  any,  in  the 
percent  reflectance  for  purities  better  than  99.95%  and, 

o 

for  wavelengths  greater  than  2200  A,  99.5%  pure 

aluminum  films  are  sufficient  to  insure  maximum  reflectance 

The  condition  of  the  oxide  layer,  which  forms 
naturally  on  aluminum  upon  exposure  to  air,  has  a  direct 
bearing  on  the  optical  properties  of  the  aluminum  coating. 
Generally,  the  oxide  layer  will  grow  (in  air)  to  a  thick- 

o 

ness  of  from  30  to  40  A  and  will  have  little  effect  on 

o 

the  reflectance  for  wavelengths  greater  than  2200  A.   At 
shorter  wavelengths,  however,  the  structure  of  the  oxide 
layer  is  most  critical,  with  a  thinner,  tighter  layer 
resulting  in  higher  reflectance  characteristics  for  ultra- 
violet radiation.   This  favorable  oxide  layer  structure 
is  obtained  by  evaporating  the  aluminum  film  at  high 

o 

deposition  rates  (not  less  than  200  A/sec)  at  the  lowest 

_5 
possible  pressure  (preferably  lower  than  10   mm  Hg) . 
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The  substrate  temperature  during  evaporation  will 
affect  the  crystal  size  (and  thus  the  surface  roughness 
and  reflectance)  of  the  film.   The  larger  the  value  of  the 
substrate  temperature,  the  lower  the  value  of  the  re- 
flectance will  be,  particularly  at  the  shorter  wavelengths 

o 

(X  <  2200  A) .   However,  if  the  substrate  temperature  is 
maintained  below  150°C  during  evaporation,  the  reflectance 
of  the  aluminum  film  will  not  be  affected  at  wavelengths 

o 

greater  than  2200  A.   A  similar  effect  results  from 
increasing  the  film  thickness  beyond  that  required  for 
maximum  reflectance  which,  for  aluminum  films,  is  a 

o 

thickness  between  700-900  A  for  wavelengths  greater 
than  1000  A. 

One  of  the  more  important  factors  influencing  the 
performance  of  thin  film  aluminum  reflectors  is  their 
age  and  the  environmental  conditions  under  which  aging 
occurs.   While  there  is  a  significant  decrease  in  the 
percent  reflectance  with  increasing  film  age  at  the 

o 

shorter  wavelengths,  such  is  not  the  case  for  X  >  2000  A. 
As  reported  by  Madden,  the  reflectance  of  aluminum  films 
(for  X  >  2200  A) ,  stored  for  longer  than  two  years  under 
normal  laboratory  conditions,  remained  unchanged.   Regard- 
less of  the  age  of  the  film,  however,  certain  environmental 
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conditions  can  cause  a  change  in  its  optical  properties. 
In  particular,  the  relative  humidity  and  exposure  of 
the  film  to  ultraviolet  radiation  can  produce  significant 
changes  in  the  reflectance  characteristics;  however,  as 
reported  by  Hass,  et  al.     ,  these  changes  are  minimal 

o 

for  wavelengths  longer  than  2000  A. 

In  reviewing  the  previous  discussion  it  is  readily 
apparent  that  the  reflectance  characteristics  of  aluminum 

o 

films,  for  X  >  2000  A,  remain  uneffected  and,  as  shown 
by  Hass,  et  al.     ,  relatively  constant  over  a  wide 
range  of  environmental  conditions.   On  this  basis,  and 
subject  to  the  provision  that  the  optical  characteristics 
of  the  film  do  not  measureably  change  upon  its  exposure 
to  the  plasma,  the  thin  film  aluminum  coating  was 
selected  for  the  reflecting  surface  of  the  radiation 
gage. 

Since  the  surface  quality  of  a  thin  film  is  dic- 
tated by  the  quality  of  the  substrate  surface,  special 
care  was  taken  in  the  preparation  of  the  cylindrical, 
copper  surfaces  prior  to  the  deposition  of  the  aluminum 
films.   Specifically,  these  surfaces  were  ground,  polished, 
and  cleaned  according  to  the  following  procedure. 

Grinding: 

i)   300  grit  silicon  carbide  paper  with 
water  lubricant. 
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ii)  600  grit  silicon  carbide  paper  with 
water  lubricant. 
Polishing: 

i)  6.0  micron  diamond  paste  (Metadi  No. 
40-6172)  with  oil  lubricant, 
ii)  0.25  micron  diamond  paste  (Metadi  No. 
40-6112)  with  oil  lubricant, 
iii)  0.05  micron  aluminum  oxide  powder  with 
water  lubricant. 
Cleaning: 

i)  ultrasonic  cleaning  with  detergent  agent 
(ALCONOX) . 
ii)  Ultrasonic  cleaning  with  isopropyl 

alcohol  -  blown  dry  with  nitrogen  gas. 

Subsequent  to  this  procedure,  the  aluminum  films 
were  deposited  with  the  vacuum  deposition  system 
(CVE-15  by  Consolidated  Vacuum  Corporation)  shown  in 
Figure  A. 2.   In  an  attempt  to  gain  a  uniform  deposit 
over  the  entire  semicylindrical  (0.5  cm  radius)  surface 
of  the  half-segment,  a  special  jig  was  designed  to 
simulate  a  line  source  of  evaporation.   As  shown  in 
Figure  A. 3,  this  device  positions  the  half-segment  such 
that  the  filament  (or  evaporation  source)  extends  along 
the  axis  of  the  segment  wall  surface.   The  filament  is 
placed  in  spring-tension  so  as  to  maintain  its  proper 
positon  upon  heating. 
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Figure  A. 2   Vacuum  Deposition  System 
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Due  to  the  close  proximity  of  the  substrate  and 
the  filament,  there  was  some  difficulty  in  controlling 
the  thicknesses  of  the  aluminum  films.   Therefore  in 

o 

order  to  insure  that  the  desired  700-900  A  film  thick- 
nesses were  achieved,  a  glass  micro-slide,  located  at 
the  same  position  from  the  filament  as  the  copper  sub- 
strate surface,  was  deposited  with  an  aluminum  film  at 
the  same  time  that  each  of  the  half -segments  were  film 
deposited.   These  glass  slides  were  then  used  in  a 
Twyman-Green  interferometer  to  determine  if  the  desired 
film  thicknesses  were  obtained.   All  of  the  aluminum 

o 

films  were  deposited  at  a  rate  of  approximately  250  A/sec 
at  a  pressure  of  2-5  x  10    mm  Hg. 

While  the  complete  spectral  distribution  (for 

o 

X  >  2000  A)  for  the  monochromatic  reflectance  (p^)  of 
the  aluminum  film  was  not  measured  in  this  study,  p^  was 
measured  at  several  wavelengths  for  comparison  with  the 
monochromatic  reflectance  data  reported  by  Drummeter  and 
HassL    .   An  aluminum  thin  film  was  deposited  on  a  water- 
cooled,  copper  disc  (shown  in  Figure  A.l)  and  this  disc  was 
employed  for  the  reflectance  measurements  of  this  study. 
These  measurements,  which  were  conducted  before  and  after 
exposure  of  the  disc  to  a  free  argon  plasma  jet,  were  made 

o 

over  the  combined  wavelength  range  (4880  to  6  32  8  A)  of  the 
argon  and  helium-neon  laser  sources  employed  for  these 
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experiments.   The  results  of  these  experiments  indicate 
that  there  is  no  measureable  change  in  the  reflectance 
of  the  aluminum  film  after  its  exposure  to  the  plasma 
environment.   At  each  of  the  wavelengths  investigated, 
the  measured  reflectance  of  the  aluminum  film  was  0.98  - 
a  value  which  is  some  6-7%  higher  than  the  corresponding 
values  measured  by  Drummeter  and  Hass.   While  the  re- 
flectance measurements  reported  by  Drummeter  and  Hass 

o 

were  averaged  (over  the  wavelength  interval,  2200  A  _< 

o 

A  _<  70,000  A)  and  used  to  determine  the  average  absorp- 

tance  of  the  reflecting  surface  (a„  =  1  -  p    =  0.07), 

2        ave 

it  is  noted  that  use  of  the  reflectance  measurements  of 
this  study  would  result  in  values  for  the  radiative  wall 
heat  flux  which  are  within  the  experimental  uncertain- 
ties shown  in  Figure  IV. 3. 

A. 5   Presentation  and  Interpretation  of 
Radiative  Wall  Heat  Flux  Measurements 

As  discussed  in  Chapter  II,  it  was  initially  felt 
that  the  presence  of  wall  contamination  in  the  arc  heating 
region  would  preclude  any  meaningful  radiation  measurements 
in  this  section  of  the  duct.   For  this  reason,  the 
radiation  gage  was  first  used  in  the  field-free  region 
(downstream  from  the  anode),  since,  in  this  section,  wall 
deposits  had  not  been  observed.   With  reference  to 
Equations  A. 7  and  A. 9  and  using  values  of  a,  =  0.98  and 
a?  =  0.07,  it  was  determined  from  the  calorimetric 
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measurements  performed  in  these  experiments  that,  of  the 
total  wall  heat  flux  for  the  reflecting  half-segment 
(q2) ,  2%  is  contributed  by  the  radiative  component. 

While  attempts  to  eliminate  the  wall  contamination 
in  the  arc  heating  region  had  failed,  radiation  experi- 
ments were  nonetheless  performed  in  this  region  in  the 
hope  that  some  meaningful  information  could  be  acquired. 
Upon  disassembly  of  the  duct  after  some  four  hours  of  high 
current  (I  >  100  amp)  operation,  it  was  observed  that, 
while  the  black  surface  of  the  radiation  gage  and  the 
copper  surfaces  of  the  adjoining  segments  were  contamina- 
ted, the  reflecting  surface  of  the  radiation  gage  was 
completely  free  of  contaminants.   Since  the  deposition 
of  contaminants  on  the  black  surface  precluded  specifi- 
cation of  its  optical  properties,  the  radiation  gage 
could  no  longer  be  used  in  conjunction  with  Equation  A. 7 
to  provide  measurement  of  the  radiative  wall  heat  flux 
in  the  arc  heating  region.   In  contrast,  the  absence 
of  contamination  on  the  reflecting  surface  continues  to 
permit  specification  of  its  optical  properties,  speci- 
fically a_.   It  was  therefore  assumed  that  the  radiation 
gage  measurements  performed  in  the  field-free  region 
could  be  used  to  delineate  the  radiative  and  convective 
wall  heat  flux  components  in  the  arc  heating  region.   That 
is,  of  the  calorimetrically  measured  total  wall  heat 
flux  for  the  reflecting  half-segment,  2%  is  due  to 
radiation. 
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On  this  basis,  a  series  of  experiments  were  con- 
ducted with  the  radiation  gage  to  determine  the  radiative 
and  convective  wall  heat  flux  for  the  asymptotic  region 
in  argon.   As  indicated  in  the  previous  section,  these 
experiments  were  further  used  to  validate  the  assumption 
that  the  presence  of  the  radiation  gage  does  not  signifi- 
cantly alter  the  local  thermal  state  of  the  plasma.   The 
results  of  this  experimental  investigation  are  presented 
in  Table  A.l. 

Prior  to  an  interpretation  of  the  data  presented 
in  Table  A.l,  it  is  noted  that  meaningful  radiation 
measurements  for  arc  currents  less  than  100  amp  could 
not  be  obtained.   This  restriction  is  due  to  the  fact 
that,  at  the  lower  arc  currents  in  argon,  the  radiative 
contribution  to  the  total  wall  heat  flux  is  too  small 
to  permit  its  resolution  within  the  estimated  uncertainty 
of  these  measurements  (±  10%) . 

The  asympotic  convective  wall  heat  flux  values 
shown  in  Table  A.l  were  obtained  by  taking  98%  of  the 
calorimetrically  measured  values  for  the  total  wall  heat 
flux  to  the  reflecting  surface.   The  first  three  of  the 
five  separate  experimental  runs  that  were  performed 
represent  an  attempt  to  determine  if  there  exist  both 
anamolies  in  the  AT  -  transducer  outputs  and  asymmetric 
conditions  within  the  arc.   While  the  data  for  the  first 
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three  runs  suggest  that  these  conditions  are  present,  the 

evidence  is  not  conclusive  since,  with  the  exception  of 

the  150  amp  value  for  Run  No.  2,  all  of  the  q       values 

^conv,a 

agree  to  within  the  experimental  uncertainty  for  wall 

heat  transfer  measurements  (±  6%)  . 

The  remaining  two  experiments  (Runs  No.  4  and  5) 

employed  a  ref lector-versus-ref lector  surface  combination 

for  the  radiation  gage,  in  which  case  the  q       results 

3  '  ^conv,a 

were  obtained  by  taking  9  8%  of  the  combined  measured 
total  wall  heat  flux  values  for  both  half-segments.   A 
comparison  between  the  data  obtained  in  these  two  experi- 
ments indicates  that  the  purity  of  the  gas  has  an  effect 
on  the  measured  results;  however,  as  was  the  case  in  the 
first  three  runs,  the  data  agree  to  within  the  experi- 
mental uncertainties. 

In  comparing  the  combined  results  for  all  five 
experimental  runs,  it  is  noted  that,  again  with  the 
exception  of  the  15  0  amp  value  for  Run  No.  2,  the  data 
agree  to  within  ±6%.   In  particular,  it  is  observed  that 
such  agreement  exists  between  the  data  acquired  with  the 
black-versus-ref lector  surface  combination  and  that 
acquired  with  the  ref lector-versus-ref lector  surface 
combination.   This  fact  suggests  that  the  presence  of  the 
radiation  gage  does  not  significantly  alter  the  local 
thermal  state  of  the  plasma.   Accordingly,  the  convective 
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wall  heat  flux  data  obtained  with  the  radiation  gage  can 

be  compared  with  the  q  measurements  obtained  with  the 

standard  copper  segments  in  the  absence  of  the  radiation 

gage.   For  this  purpose,  the  q       values  for  each 
3  r  v     v        '      ^conv,a 

arc  current  condition  were  algebraically  averaged  and 
the  resulting  values  were  subtracted  from  the  corres- 
ponding q  values  (Table  IV. 1)  to  give  the  radiative  wall 
heat  flux  results  shown  in  Table  A.l.   The  values  for 
the  fraction  of  radiative  to  total  wall  heat  flux  shown 
in  Table  A.l  were  used  for  the  comparisons  in  Chapters 
III  and  IV. 
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